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HURRICANE EFFECT ON SEA LEVEL AT CHARLESTON 


Bernard D. Zetler* 
(Proc. Paper 1330) 


ABSTRACT 


_ Tropical cyclones from 1922 to 1955 have been classified according to the 
cientation of their tracks with respect to Charleston Harbor. Theoretically, 
1e displacement of sea level depends on this orientation; the tabulated dis- 
‘acements caused by tropical cyclones during 34 years of tide observations 
ear out the relationships. 


. 
. 


_ When it was decided in the Coast and Geodetic Survey to make a study of 
ie effect of hurricanes on sea level, Charleston was chosen as the area of 
ie pilot study for several reasons. Its position makes it vulnerable to peri- 
ic hurricanes, and the low elevation of its downtown section creates the 
eed for accurate forecasting of storm heights. Furthermore the problem 
i a to be simpler for a port just off the ocean than, for example, a port 
a Chesapeake Bay in which interference oscillations may be anticipated. 
astly there was available for Charleston a relatively long series of tide 
»servations dating back to 1922. 
The attempt in this study has been to classify hurricanes into categories 
ich that some groups are likely to have a large effect on sea level while 
hers have little or none. Figure 1 shows the basis of one set of criteria 
nich was chosen. Considering the Atlantic coastline near Charleston to lie 
a SW to NE direction, the optimum direction of water motion to raise sea 
vel is from the SE. Considering the Coriolis effect in the northern hemi- 
here of water moving about 45 degrees to the right of the wind, water move- 
ent from the SE will be set up by easterly winds. If an observer faces into 
e wind, the storm center will be on a line about ten points (112°) to the 
ght. Therefore with easterly winds at Charleston, the storm center will 
ve a bearing of approximately SSW. Accordingly, all other things being 


a ee nf : 

»te: Discussion open until January 1, 1958. Paper 1330 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 1, August, 1957. 


*Oceanographer, Coast and Geodetic Survey, U. S. Dept. of Commerce, 
Washington, D. C. 
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equal, storms having a bearing from Charleston of about SSW will cause a 
rise in sea level at Charleston while storms having a bearing of about WNW 
or ESE will cause currents to move parallel to the coast with little effect on 
sea level. A tropical cyclone bearing NNE from Charleston should lower 
sea level by moving water away from the coast. Storms having other bear- 
ings from Charleston will have effects proportionate to their departures 
from the above directions. 

Figure 2 shows eight octants radiating out from Charleston, each repre- 
senting an arc used in classifying the bearing of a cyclone center. Each 
name represents the direction in the center of the octant. By tabulating the 
octants through which the center of a cyclone moves in approaching and mov- 
ing away from Charleston, one really is enumerating the changing bearing of 
the storm center as observed from Charleston. 

Those storms classified as having Charleston on the right side of their 
tracks ordinarily move inland whereas those with Charleston on the left side 
remain out at sea. The criteria was first considered important on the as- 
sumption that storms moving inland pass a critical depth in shallow water at 
which resonance contributes to a surge.(1) Using this reasoning in studying 
the surge effect at Charleston, it would be necessary to distinguish between 
storms crossing the Atlantic coast and those crossing the Gulf coast. How- 
ever this study shows large positive displacements with the track crossing 
the coastline in the Gulf of Mexico as well as the Atlantic coast. Apparently 
the criteria of Charleston being on the right side of the track is important 
because of associated strong winds over shallow water outside the harbor 
rather than because of a surge being set up at a critical depth. A storm with 
Charleston on the left side of its track contributes to changing sea level in 
the form of forerunners which are moderate rises in sea level over a period 
of time due largely to currents set up by wind circulation. In deep water 
much of this effect is dissipated before reaching the coast. 

The question of how close a cyclone had to come to Charleston to affect 
sea level appreciably was resolved by including all storms shown on the 
track charts for the first few years considered. It became apparent that it 
was possible to restrict the study to cyclones coming within 345 miles (300 
nautical miles) of Charleston. This is a convenient figure to work with as it 
is measured on a chart by five degrees of latitude near the latitude of 
Charleston. In the final tabulations, only the portions of tropical cyclones 
within a radium of 345 miles of Charleston are considered. This probably — 
includes some cyclones which need not be considered but the effective 
diameters of tropical cyclones vary considerably in size and the diameters © 
of particular storms are not usually listed in the published statistics. 

The cyclones are classified first according to the octants in which they — 
came closest to Charleston but the data on other octants are also important. 
For example two cyclones may have been closest to Charleston in Octant ES 
If one of them has passed through Octants SSW and SSE in approaching ESE, 
it presumably is more likely to have raised sea level than the other. The _ 
final tabulation also subdivides the storms classified within an octant to ~ 
those with Charleston on the right side of the cyclone track and those with — 
Charleston on the left. Tropical cyclones with Charleston on the right side 
are presumed to have a greater effect on sea level. The results bring this 
out very clearly. . 
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Figure 1. Theoretical bearing of center of tropical cyclone for maximum 
rise in sea level at Charleston. 
. 


Figure 2. Octants radiating from Charleston used in classifying 
the bearing of the center of a tropical cyclone. 
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Tide Records 


In determining the displacement of water level caused by a tropical cy- 
clone, it was found desirable and necessary to consider the continuous record 
as represented by the complete tide curve. In most cases, a comparison of 
the observed and predicted heights on each hour presented a true picture, but 
it was found unsatisfactory for final results to compare heights only at the 
times of high and low tides, for, being so far apart in time, they seldom 
showed the maximum variation. Furthermore the time of high or low tide 
may be advanced or retarded several hours by a storm and a comparison of 
observed and predicted high or low tides would deal with water level at dif- 
ferent times. 

Despite these limitations comparisons of observed and predicted high and 
low tides were used as the preliminary step in determining what storms were 
to be included in the study. Any storms for which the calculated displace- 
ments did not include at least one value of one foot or more were eliminated 
from further consideration. These are shown in the tabular results as having 
a displacement of less than one foot. Obviously, in using this arbitrary rejec 
tion limit, some storms in which the calculated maximum hourly displace- 
ment would have been somewhat over a foot were rejected from further study 
However, in each case the observed tide curve was examined for unusual 
fluctuations between the high and low waters in order that no important dis- 
placements would be rejected by this method. 

Predicted tides at Charleston are referred to a datum of mean low water, 
based on gage records for the 19-year period 1924 - 1942, and take into ac- 
count seasonal variations in sea level. Inasmuch as sea level along our east 
coast has been rising gradually (approximately 0.4 foot in the last few dec- 
ades), each year’s observations were referred to mean low water observed 
in the same year before being compared with predictions. 


Tabulation of Tropical Cyclone Data and Sea Level Displacement 


The following tables furnish data on tropical cyclones from 1922 to 1955, 
describing their paths, speeds, and their bearings and distances from 
Charleston. They indicate the maximum displacement of sea level and the 
interval of time between the cyclone coming closest to Charleston and the 
maximum displacement. When the maximum displacement occurred before — 
the cyclone reached its closest point, the time lag is preceded by a negative - 
sign. The maximum water height (referred to Mean Low Water) attained 
during the period of the storm is given in the last column. Underlined names 
of storms indicate that they were hurricanes. 

Tropical cyclones are numbered in accordance with Weather Bureau trail 
charts except that those earlier than 1933 are numbered in accordance with | 
the track charts in “HURRICANES” by I. R. Tannehill. Numerous adjectives 
were found to describe the intensity of tropical cyclones in various sources. 
Even the clear-cut distinction of whether or not a tropical cyclone attained - 
hurricane status is somewhat unsatisfactory since that degree of intensity 
may not apply to the storm during the period it may have affected sea level — 
at Charleston. Some track charts show the position of the storm center ; 
every twelve hours; others give it every twenty-four hours. It was necessar- 
to infer that the center of the storm moved at a uniform rate between adjace 
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Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant SSE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum 
Length Distance Of Sea Level 
Tropical of to Time 
lone Octant Fath Duration Speed Charleston Height La 
miles hours m.peh. miles feet hours 


| Charleston on right side of cyclone track 


20 II 

SSW 
1951 Able ESE 
1st Approach) SSE 


Charleston on left side of cyclone track 


| 
938 vII WwSW 


SSE 110 11 10 
235 18 13 
215 li 20 
220 na 20 


205 4 51 

) ssw 150 3 50 
) SSE 20, (O04), (50) 
ESE 20 (0.4) (50) 

ENE 340 (7) = (49) 

24 IV WSW =-185 19 10 
ssw 145 9 16 

SSE 40 3 13 

ESE 55 2 28 
| ENE 275 i 25 
1932 Iv SSW -285 1 26 
SSE 65 3 22 

ESE 95 4 2k 

| ENE 270 lk 23 
p37 I yee ees ea 
SSE 105 6 18 

| ESE 105 6 18 
ENE 265 15 15 
953 Hazel SSW 110 5 22 
| SSE 195 7 28 
| ESE 215 8 27 
38 VI SSE 140 ? ? 

sl III SSE 275 3h 8 
st Approach) ESE 45 3 15 
31 IV ssw 60 6 10 
SSE 195 21 a 
ESE 100 ? ? 

SSE 195 1 18 

ESE 100 4 25 


140 
170 


235 
230 


+2,.2 


+1,.1 


<1.0 


<1.0 


<1.0 


<1,0 


+123 


<1.0 
<1.0 


<1.0 


<1.0 


=22 
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Displacement Maximum 


Height 
Above 
MLW 

feet 


Tek, 


509 


6.8 


6.6 


beh 


505 


6.7 


6.9 
6.1 


6.2 


508 
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: 
TABLE 2 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant SSW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimun Displacement Maximur 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height lag MLW 
miles hours mepeh. miles feet hours ee 
Charleston on right side of cyclone track 


1940 III ESE 100 13 8 280 +5.1* 2? - 200% 
SSE 250 32 8 35 
SSW 30 5 6 30 
WSW os 15 6 35 
WNW 250 23 Gh 125 
1947 IX SSE 5 6 16 re 42.6 4 8.5 
ESE - 300 
) Course (46 ~~ (215) 
ESE ) unce rtain, a 110 
SSE 75 8 9 80 
ssw 70 6 12 70 
WSW 95 iz - 100 
1947 VIII SSE 210 (13) (16) 16 +1.3 2 6.7 
SSW 150 7 21 140 
WOW 200 ? = 160 
1936 IX SSE 125 6 21 270 <1,0 - 567, 
SSW 215 12 18 24,0 
ww 125 (9) (14) 275 
1937 III SSE 85 & 1 300 <1,.0 - 6.6 
SSW 230 22 10 250 
WSW 45 wy 10 275 
1939 II SSW 210 16 13 315 <1.0 - 6.7 
(1st Approach) WSW 40 5 8 330 


Charleston on left side of cyclone track 


1946 VI SSW 285 ? ? 60 +1,4 ? 602 
1950_XII SSW 140 ? - 230 +167 ? 607 
1937 VI WwSW 90 12 8 26 <1.0 = 

poe! ps & a 1. 6.5 
1941 VI SSW 220 ? - 240 0 - 7.0 


* According to measured high water mark, At time gage ceased ; 
displacement was +3.7 feet, B88 operating, 
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TABLE 3 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant WSW When Closest to Charleston 


: 
' 
. 
. 
Storm Data for Octants Traversed 


Maximum 
Minimun Displacement Maximum 
Length Distance of Sea Level Height 
Tropical of to Time Above 
Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 
miles hours m.p.h. miles feet hours feet 
Charleston on right side of cyclone track 
1953 Florence WSW 285 ? ? 40 +105 -72 Tel 
1934 I SSW 70 4 18 300 +309 =3 9.1 
SSE 145 9 16 130 
SSW 90 8 2st 45 
. WSW 40 3 13 40 
WNW 45 6 8 50 
NNW 80 7, 5 115 
WNW 150 % = 95 
1927 V SSE 145 8 18 205 (41.9 +1 5.3 
| ssw 150 10 15 aol {20 +8 
WSW 50 3 17 60 
WNW 65 4 16 65 
NNW 225 (12) (19) 85 
NNE 65 ( 4) (16) 270 
| 1952 Able SSE 205 18 ila 190 +3,0% +2 6.6 
ssw 155 12 13 70 
) vSW 55 3 18 60 
| WNW 70 4 18 65 
NNW 215 16 13 105 
NNE 80 ( 6) (13) 270 
1950 V SSW 330 60 6 180 +2.7 +5 6.8 
| ww 190 13 15 175 
| WNW 85 6 14 270 
1926 I ssw 255 13 20 215 +2.8 +42 77 
po wow 195 10 20 210 
WNW 60 5) 20 265 
1950_XI ssw 230 18 13 230 +23 +6 609 
: ‘ WsW 230 ? - 225 
1933 VI ssw 45 6 8 320 <1,0 - 6.6 
WSW 180 2 = 235 
1953 Alice WSW 85 ? ? 300 <1,.0 - 509 
1948 III WSw 60 11 5 330 <1.0 - 703 
: WNW 45 7 6 335 


- The maximum displacement may have been exceeded and the time lag changed 
during the period the gage was inoperative. 
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TABLE ha 
Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant WNW When Closest to Charleston 
Storm Data for Octants Traversed 

Maximum 
Minimun Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 
miles hours mph. miles feet hours feet 


Charleston on right side of cyclone track 


1924 V WSW 320 1 29 30 (4259  o8 6.7 
WNW 25 al 25 25 (-2.2 +5 
NNW 25 A 25 25 
NNE 310 (10) (31) 35 

1945 IX SSW 260 26 10 110 +10 3 709 
WSW fe: 5 15 60 
WNW 50 3 17 50 
NNW 90 6 15 60 
NNE 205 12 aly 135 

1944 X SSW 255 17 its 105 +3.8 =3 8.5 
WSW 85 4 21 65 
WNW 45 2 22 60 
NNW 100 5 2 65 
NNE 240 ual 22 130 

1947 VII SSW 195 12 16 180 +262 a) 605 
WSW 130 7 19 75 
WNW 45 2 22 60 
NNW 60 3 20 65 
NNE 275 13 al 95 

1946 V SSW 215 10 22 150 42.3 =, 702 
WSW 105 5 al 90 
WNW 16 3 25 85 
NNW 100 8 12 90 
NNE 205 15 125 
ENE 100 ( 6) 17 270 

203550 SSW 
WSW 
WNW 90 - - a 
NNW 
NNE 

The tide gage at Charleston was inoperative during this storm. 

1933 XII SSW 200 28 7 270 +10 (-40)* (6,3). 
WSW = 230 2h 10 125 ~ 
WNw 80 ? - 100 . 
NNW 40 2 - 105 


* The maximum displacement may have been exceeded and the time chan ed 
during the reriod the gage was inoperative. a 
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TABLE 4b 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant WNW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 


miles hours mp.h. miles feet hours feet 


Charleston on ht side of cyclone track = cont'd 


| 1949 II SSW 185 uy, 13 24,0 Sy ea a 765 
| WSW 165 1 15 150 (-1.5 +10 

WNW 110 5 22 135 
| NNW 230 9 26 150 
| NNE 50 1 50 305 
1928 IT WSW 235 20 12 310 <1.0 - bo4 
ww = 260 21 12 300 

NNW 100 4 uy 305 


1932 VII WNW 310 ~ - = 
NNW 
The tide gage at Charleston was inoperative during this storm. 
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TABLE 5 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant NNW When Closest to Charleston 


Storm Data for Octants Traversed 


Maximun 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Fath Duration Speed Charleston Height Lag MLW 
miles hours m.peh. miles feet hours feet 


Charleston on right side of cyclone track 


1928 V SSW 295 19 16 90 (43.8 -6 (8.1) 
WSW 90 6 15 30 (-1.2 +21 
WNW 20 2 10 20 
NNW 20 2 10 15 
NNE 330 2k lu, 20 
1941 V WwSsW 320 26 12 70 <1.0 - 57 
WNW 70 5 uy, 35 
NNW 25 2 12 20 
NNE 30 2 15 25 
ENE 5 9 13 25 
ESE 235 18 13 130 
1929 II Wwsw = 30 24 13 90 (+1.6 +1, 7.0 
WNW 70 5 uy 65 (-2.2 +10 
NNW 80 5 16 60 
NNE 250 12 21 100 
1928 I ssw 220 32 7 285 +14 6 6.3 
WSW = 230 2h 10 225 
WNW so155 12 13 155 
NNW 140 10 uy 145 
NNE 200 10 20 190 
1936 XVI wsw 170 n 15 <1.0 - 6.3 
www =s:1170 uy 12 165 
NNW  =160 6 27 160 
NNE 160 6 27 220 
1939 II WNW 200 18 ala} 270 <1.0 - 567 
(2nd Approach) NNW 230 12 19 250 
NNE 35 13 23 310 
1934 II NNW =: 120 6 20 310 <1.0 = — 5.39 
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TABLE 6 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant NNE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Lag MLW 


miles hours mepeh. miles feet hours feet 


Charleston on left side of cyclone track 


1942 VI ENE 50 4 12 310 <1.0 - 70 
| NNE 115 ? - 2h0 

: 
| 
points to obtain a time lag measured in hours. Occasionally an extrapolation 
| was necessary for this measurement and parentheses are set around these 
figures to indicate an inference. A number of cyclones dissipated near the 
points on their tracks closest to Charleston; a time lag estimate would be un- 
reliable for these and a question mark appears in the column. Within each 
classification, the tropical cyclones are listed in order of increasing mini- 
mum distance to Charleston. Several cyclones came within the 345 mile 
range of Charleston, moved away, and then came back again. These ap- 
proaches are listed independently in the tables but are distinguished by a 
Rotation “1st approach” or “2nd approach.” 


Tropical Cyclones In Octant SSE 


_ Theoretically a cyclone in Octant SSE will raise sea level but not as much 
as one in SSW. However a cyclone approaching Charleston from SSW will 
usually pass over Florida and in doing so, lose some of its energy. 

- There are 11 storms listed in this classification. There was a rise in 
sea level for the two cyclones with Charleston on the right side of their 
cyclone tracks. All nine having Charleston on the left showed little dis- 
placement. 

Hurricane II, 1940, had a displacement of +2.2 feet 24 hours after coming 
closest to Charleston. This storm moved from Octant SSE into Octant SSW 
and the maximum displacement took place while the storm was in SSW, ac- 
counting for the large time lag. 

Hurricane Able, 1951, had a maximum displacement of +1.1 feet, two hours 
before coming closest to Charleston in its first approach. It then turned 
counter -clockwise, out of range of Charleston, without crossing the coast, 
thus accounting for the small displacement and small lag. Hurricane II, 1940, 
Came almost a hundred miles closer to Charleston than did Hurricane Able, 
1951, and moved considerably slower. 

_ Four of the remaining nine came within less than one hundred miles of 
Charleston and all four passed through Octant SSW. One of them, IV 1924, 
attained hurricane intensity but had a small diameter. The fact that 


k 
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TABLE 7 


Storm Tides at Charleston, Se Ce, 1922 to 1955 
Tropical Cyclones in Octant ENE When Closest to Charleston 


Storm Data for Octants Traversed 


Maximum 
Minimum Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Ia MLW 
miles hours mepeh, miles feet hours feet 


Charleston on left side of cyclone track 


1954 Hazel SSE 220 8 28 145 (+127 = Tel 
ESE 105 rm 26 100 (1.5 15 
ENE 105 3 35 90 
NNE 255 5 51 125 

1944 IIT SSE 165 ln 15 180 <1.0 ~ 6.5 
ESE 140 9 16 105 
ENE 80 5 16 100 
NNE 230 13 18 120 

1 Connie ESE 255 40 6 165 +2,.0 -7 609 
ENE 200 32 6 160 
NNE 60 4 15 315 

1955 Diane ESE 215 1, 15 185 ped ee +5 7.8 
ENE 145 12 12 170 
NNE 230 17 14 185 

1953 Barbara SSE 30 4 8 290 <1,0 - 5.8 
ESE 230 25 9 185 
ENE 195 15 13 180 
NNE 35 «! 12 265 

1955 Ione ESE 230 13 1s 220 +1,.3 -10 6.9 
ENE 160 13 12 215 
NNE 90 10 9 255 

1944 VI ESE 230 14 16 230 <1.0 - bel 
ENE 235 10 2h, 225 

1 XIII ESE 140 15 9 270 <1.0 - 509 
ENE 195 18 unk 265 

1936 XIII ENE 185 16 12 310 <1.0 - 6.18 

1940 IV ESE 35 4 9 330 <1.0 - 6.9 - 
ENE 145 15 10 310 3 

1933 VIII ENE 105 8 13 325 41.7 im2uie Gehl 
NNE 50 4 12 330 

1 III ESE 50 6 8 330 41.8 =—2 Were 

2nd Approach) ENE 100 1 9 325 


ASCE ZETLER 1330-13 


Charleston was on the left side of their hurricane tracks appears to be the 
explanation for their failure to raise sea level appreciably. 


| Tropical Cyclones in Octant SSW 


Theoretically, cyclones with centers in Octant SSW will have the maximum 
ability to raise sea level at Charleston. Of the six cyclones which had 
_ Charleston on the right side of their tracks, the three closest caused notable 
rises in sea level; the three furthest away did not. Furthermore the amount 
of displacement by the first three appears directly related to their minimum 
distances to Charleston. This.ordered sequence may also be due to the first 
two cyclones, III 1940 and IX 1947, attaining hurricane intensity while the 
third, VIII 1947, did not. The time lags for the second and third cyclones are 
negligible, as is to be expected of most cyclones coming closest to Charleston 
in Octant SSW. The time lag on the first, with a displacement of +5.1 feet, 
_ was not determined because of tide gage difficulties. 
| The loss of some tide data for Hurricane III 1940 illustrates a basic diffi- 
_Culty in a study such as this. The storms which strike Charleston hard and 
Cause large displacements in sea level frequently cause tide gage failure and 
On a few occasions have destroyed the tide gage installation. The loss of 
| record of some extreme conditions is a serious handicap to a study such as 
this. The program now being carried out of modifying the tide gage design, 
_Taising the gages and remoting recorders will help in the future. 
| There were four cyclones which were closest to Charleston in Octant SSW 
: with Charleston to the left of their tracks. Tropical cyclone VI 1946, which 
| 


Came within 60 miles of Charleston caused a displacement of +1.4 feet. Hur- 
ricane XII 1950, was 230 miles from Charleston and caused a displacement of 
+1.7 feet. The other two, both 240 miles from Charleston, had little effect on 

sea level at Charleston. Since all four dissipated in Octant SSW there is no 

data on time lag. 


Tropical Cyclones in Octant WSW 


As is to be expected in view of the usual paths of West Indian tropical cy- 
clones, all the storms which came closest to Charleston in Octant WSW have 
Charleston on the right side of their cyclone tracks. Since many of these 
passed through Octant SSW and even SSE prior to reaching WSW, a large 
positive sea level displacement with little time lag is to be expected from 
these cyclones. This is found to be the case. 

Six of the ten cyclones in this category raised sea level from 2 to 4 feet; 
the time lag in each case was small. Of the remaining four, one dissipated in 
the octant and the other three were not hurricanes and their distances from 
Charleston were large. 

Sometimes the elevation of sea level, as caused by a tropical cyclone, is 
followed by a lowering of sea level appreciably below the predicted value. 
The maximum negative displacement sometimes occurs just a few hours after 
the maximum positive displacement; at other times the time difference is a 
day or more. As noted previously, tropical cyclones in Octant NNE should 
lower sea level at Charleston most effectively and cyclones in Octants NNW 
and ENE should have similar effects. None of the tropical cyclones classi- 
fied as coming closest to Charleston in Octants SSE and SSW passed through 
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TABLE 8a 


Storm Tides at Charleston, S. C., 1922 to 1955 
Tropical Cyclones in Octant ESE When Closest to Charleston 


Storm Data far Octants Traversed 


Maximum 
Minimim Displacement Maximum 
Length Distance Of Sea level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Ia MLW 
les hours mepeh. miles feet hours feet 


Charleston on right side of cyclone track 


1934, III ENE 120 ? - 60 “1.2 -16 5.6 
ESE 50 8 6 50 
SSE 140 16 9 60 
SSW 255 17 15 155 

1935 V ENE 100 10 10 265 +1e4 +413 bok 
ESE 230 23 10 240 
SSE 70 7 10 300 
Charleston on left side of cyclone track 

1930 II SSW 230 LA 5 160 <1,0 - 5.8 
SSE 115 12 10 10 
ESE 115 12 10 100 
ENE 235 18 13 145 

1945 Tt SSW 265 21 13 140 42.6 +1 7.0 
SSE 100 6 ale: 105 
ESE 100 6 17 100 
ENE 265 4, 19 125 

1925 III SSW 160 6 27 230 +2,.0 +2 609 
SSE 165 aah 15 120 
ESE 110 6 18 110 
ENE 165 13 13 130 
NNE 110 9 12 230 
ESE 260 72 if 140 
ENE 255 15 17 160 

1934 VI SSE 190 L7. nu 185 <1.0 - 5.9 
ESE 1,0 6 23 150 
ENE 230 8 29 170 

1922 III ssw 60 ? - 290 +1,8 + Tee. 
Sse 195 (29) (7) 185 “7 
ESE 190 39 5 180 
ENE 100 18 6 260 

1951 How SSW 50 3 17 320 +1,.2 nfl Te2 
SSE 215 12 18 215 . 
ESE 180 12 15 195 | 
ENE 150 18 8 270 
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TABLE 8b 


Storm Tides at Charleston, S. Ce, 1922 to 1955 
Tropical Cyclones in Octant ESE When Closest to Charleston 


| Storm Data for Octants Traversed 


Maximum 
. Mininun Displacement Maximum 
Length Distance Of Sea Level Height 
Tropical of to Time Above 


Cyclone Octant Path Duration Speed Charleston Height Lag MLW 


miles hours mepeh. miles feet hours feet 


| Charleston on left side of cyclone track - cont'd 
| 


1933 XVI SSE 15 Or ey 330 41.3 (+14) 609 
| ESE 205 ? - 200 
1937 II SSE 10 13 8 240 <1.0 - 6.7 
| ESE 190 18 uk 200 
| ENE 160 12 13 220 
1924 II SSE 195 40 5 265 +165 ~-7 6.1 
| ESE 195 12 16 210 
i ENE 185 8 23 235 
1950 I ESE 235 25 9 225 <1.0 = bok 
: ENE 180 10 1g 230 | 
1948 IV ESE 180 19 9 21,0 <1,0 = 70 
: ENE 125 Boe a. 10 245 
1949 I ESE 230 15 15 240 <1.0 =“. ° 6,8 
: ENE 200 10 20 245 
: ENE 15 v) a1 255 
1954 Edna ESE 260 28 9 250 G13 +2 609 
7 ENE 175 10 18 265 (41.3 +15 
1937 VIII SSE 130 9 wy 275 <1.0 = 602 
ESE 230 11 a 250 
ENE 80 3 27 280 
938 9 27 255 Tide gage not oper= 
<8 oe 185 5 37 260 ating at this time, 
1944 1 SSE 40 4 10 325 <1,0 - 607 
ESE 310 18 6 275 
1951 Able ESE 255 21 12 275 <1.0 - Toh 
(@d Approach) ENE 160 12 13 280 
1943 II ESE 200 14 an 290 <1.0 - Teh 
: ENE 65 5 13 310 
1932 IIT SSE 115 16 7 305 <1.0 - 6.5 
ESE 195 20 10 300 
1951 Jig ESE 160 15 nl 325 41,5 =15 700 
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Octants NNE or NNW, a few of them were in ENE but far off from Charleston. 
Consequently none of them lowered sea level at Charleston appreciably. 
Tropical cyclone V 1927 caused a negative displacement of sea level of 
1.0 foot seven hours after a positive displacement of 1.9 feet. The negative 
displacement occurred while the cyclone center was in Octant NNW although 
it continued into NNE; however, the minimum distance to Charleston in NNE 
was 270 miles so the storm may have been out of effective range of Charles- 
ton by the time it reached NNE. Two other cyclones I 1934 and Able 1952 
passed through NNW also. Neither lowered sea level by a foot but V 1927 
was closer to Charleston while in Octant NNW than either of them. 


Tropical Cyclones in Octant WNW 


As was true of tropical cyclones in the previous category, tropical cyclone: 
coming closest to Charleston in Octant WNW will ordinarily have Charleston 
on the right side of their track and therefore will cause positive displace- 
ments of sea level at Charleston because they have passed through Octant 
WSW and possibly SSW. However since they contribute little to raising sea 
level while in WNW, a negative time lag is to be expected. 

This is generally true of the cyclones listed in this category in the table. 
Four cyclones caused positive displacements of 2-1/4 to 4 feet with negative 
time lags of a few hours. Two others however, 1924 V and 1949 II, caused 
positive displacements of 2 feet with small positive time lags; these two 
were also the only cyclones in the group to cause subsequent negative sea 
level displacements. While the positive time lags may possibly be attributed 
to error introduced by inferring a constant speed between points of refer- 
ence, there is no apparent reason why only these two cyclones should cause 
negative displacements. Tropical cyclone V 1924 was very close to Charles- 
ton while in Octants NNW and NNE, but a number of the other storms were 
closer than Hurricane II 1949 while in these octants. Perhaps the unusual 
severity of Hurricane II 1949 accounts for its effect. 

The effect of Tropical Cyclone V 1924 was unusual too in that there was 
only a two hour difference in time between the maximum positive and nega- 
tive displacements of sea level. This is shown on the marigram as a re- 
versal of the tidal rise between 5 and 7 a.m. (75° W. time) on September 30. 
Immediately thereafter, the tide rose again at a somewhat accelerated rate. 
The eye of the cyclone almost passed over Charleston with strong winds first 
building up the level of the sea, the calm then releasing the accumulation and 
strong winds shortly thereafter driving the water from the coast. There was 
a change in sea level of 4.1 feet which took place in two hours that was direct 
ly attributable to meteorological conditions. 

Hurricane XII 1933 caused a positive displacement of 1.0 foot with a time 
lag of -40 hours. The unusual time lag may be due in part to tide gage failure 
during the next eleven hours during which time the maximum recorded dis- 
placement may have been exceeded and in part to weakening of the hurricane 
as it approached Charleston; the storm dissipated shortly thereafter. 

One hurricane 300 miles from Charleston had little effect on sea level. 
The Charleston tide gage was inoperative during two other storms. 
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Tropical Cyclones in Octant NNW 


Those cyclones coming closest to Charleston in Octant NNW have Charles- 

| ton on the right side of their tracks. They may have large positive displace- 

| ments of sea level with negative time lags based on their effect in previous 
octants. They may also have negative displacements with positive time lags 
based on subsequent effects in Octant NNE. 

/ Hurricane V 1928 caused a positive displacement of 3.8 feet followed bya 

_hegative displacement of 1.2 feet; the time lags corresponded to the above 
estimates. Hurricane I 1928 caused a displacement of +1.4 feet with a small 

| negative time lag. 

) Hurricane II 1929 caused a positive displacement followed by a negative 

| displacement but the time lag for the rise in sea level is +4 hours. Here 

again the apparent discrepancy may be due to inference of a uniform speed 

between two points on the travel path twenty-four hours apart. The time lag 

) for the negative displacement puts the cyclone center in Octant NNE at the 

| time. 

) Four other listed cyclones in this category had no appreciable effect on 

| sea level at Charleston. Hurricane V 1941 is among these although it came 

within 20 miles of Charleston. Its failure to raise sea level is explained by 

its wide swing around Charleston in its clockwise motion but its failure to 

subsequently lower sea level while coming close to Charleston in Octants 

| NNW, and NNE and ENE is not readily explained unless the hurricane lost its 

| strength prior to approaching Charleston. 


Tropical Cyclones in Octant NNE 


Only one cyclone IV 1942, was closest to Charleston in Octant NNE. It 
'was 240 miles away and had little effect on sea level at Charleston. 


Tropical Cyclones in Octant ENE 


| Twelve tropical cyclones, all of hurricane intensity, were closest to 
' Charleston in Octant ENE. Charleston was on the left side of the track of 
} each, 
) There was a positive displacement of between one and two feet by six of 
these hurricanes, each with a negative time lag indicating that they were set 
be while the cyclone was in the ESE direction. The negative displacement of 
(15 feet for 1954 Hazel is interesting as it illustrates the behavior of a re- 
surgence as defined by Redfield and Miller. There was a positive displace - 
ment followed consecutively by a negative, a positive and then a negative dis- 
placement, each about six hours apart. These oscillations are believed due 
to the free motion of the water in restoring the disturbed levels, augmented 
in some cases by the changed direction of the wind. The period of oscillation 
following this storm, roughly that of the semidiurnal tidal cycle, appears to 
be the usual pattern at Charleston. If this is so the resurgence at Charleston 
is not apt to attain substantially greater heights than the initial rise as it will 
occur at about the same tidal phase. 
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Tropical Cyclones in Octant ESE 


There are twenty-four tropical cyclones listed as coming closest to 
Charleston in Octant ESE. Of these only two have Charleston on the right 
side of their tracks. 

Although one of these raised sea level while the other lowered it, the data 
are not contradictory. Tropical cyclone III 1934 lowered sea level at Charles 
ton by 1.2 feet. The negative time lag indicates it caused this while in Octant 
ENE. Hurricane V 1935 raised sea level 1.4 feet. The positive time lag in- 
dicates this occurred after the storm reached Octant SSE. 

Theoretically, tropical cyclones in Octant ESE will set up currents paral- 
led to the coast and have little effect on sea level at Charleston. This may be 
particularly true of cyclones with Charleston on the left side of their tracks. 
Some of the tropical cyclones listed in this classification caused positive dis- 
placements of sea level. In general these displacements are presumed to be 
due to their effect in octants previously traversed. The negative time lag or 
small positive time lag supports this suggestion. 

The displacements attributed to several tropical cyclones are probably due 
to extraneous conditions. Sea level rose 1.3 feet when tropical cyclone XVI 
1933 was near Charleston but this maximum displacement occurred about a 
half day after the cyclone had dissipated at its closest point to Charleston. 
There was a rise in sea level of 1.5 feet near the time Hurricane Jig 1951 
was first reported but equal displacements of sea level were reported on the 
previous two days. In both cases, the displacements appear to be coincidental 
to the presence of the tropical cyclones or at best they are only partially 
responsible. 

The two storms Hurricane II 1930 and Hurricane I 1945 illustrate the dif- 
ficulty in the classification of storm intensity. Both cyclones are listed as 
attaining hurricane intensity and traveled over almost identical paths while in 
the vicinity of Charleston. The first had little effect on sea level at Charles- 
ton and the second raised sea level 2.6 feet. Further investigation disclosed 
that the first had lost most of its strength before reaching the Atlantic and the 
second retained its hurricane strength, thus explaining the difference in effect 


SUMMARY 


While it has been demonstrated that the bearing of a tropical cyclone and > 
the direction in which the cyclone is moving with respect to Charleston are 
important factors in determining sea level displacement at Charleston, too 
many other variables are involved to permit accurate quantitative forecasts 
at present. It is possible to predict that certain storms will have little or no 
effect. It is also possible to predict that other storms will raise sea level by 
an estimated amount with the extreme coming within a few hours of a certain 
time, provided, of course that there is an accurate forecast of the cyclone 
track. 
As indicated previously, the intensity of a tropical cyclone at the time it — 
can affect sea level at Charleston has not been adequately determined for the 
storms in this study. Similarly the dimensions of the effective diameters of 
most storms were not available although some data are published(2) and 
more information of this kind will result from the current National Hurricane 
Research Project. The speed of the cyclone centers in the various octants 
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are given in the table but no apparent relationship to sea level displacement 
tas been noted. A tropical cyclone moving slowly should have a greater effect 
on sea level because winds will blow from the same direction for a longer 
period of time. On the other hand, a tropical cyclone moving with rapid 
velocity (forty to fifty miles per hour) will augment the velocity of winds on 
the right side of the track and thus may strongly affect sea level. Presuma- 
bly, heavy rainfall over the adjacent land area can raise the level of water in 
Charleston Harbor temporarily. 

All these factors really need to be considered to a certain extent in fore- 
casting sea level displacement. On the other hand, a study which classifies 
each cyclone according to each criteria may find that, with a limited sample, 
each storm is in a class by itself. In view of these variables, the fairly con- 
sistent results within each classification used in this study indicate the rela- 
tive importance of the criteria used. 
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SYSTEMATIC CHANGES IN THE BEDS OF ALLUVIAL RIVERS@ 


Walter C. Carey,! M. ASCE and M. Dean Keller,2 J.M. ASCE 
(Proc. Paper 1331) 


| ABSTRACT 

. 

_ This paper describes the sand wave formations present in alluvial rivers 

and discusses their variations with respect to river stage. The related rise 

sand fall of crossings is also discussed. The illustrations are prepared from 
\fathometer surveys and show the changes in the sand waves due to changes in 
Stage. 


INTRODUCTORY 


Close association with various aspects of planning, constructing and oper- 
ating works for the improvement of the Mississippi River, has gradually led 
ig the conclusion that systematic changes take place in the river bed, and that 
hese changes are a major factor in determining resistance to flow and in 
pene the general behavior of the river. Two somewhat over-lapping 


) 
| 


henomena are involved; namely, (1) a rise and fall of the crossings between 
nds in sympathy with variations in discharge, and (2) a system of sand 
waves on the river bed, which varies as to general arrangement, and as to 
size and spacing of individual waves, in accordance with some controlling 
condition. 
As to the first phenomenon, the variation in height of crossings with re- 
| to discharge, it will be noted that its existence is well known to practi- 
al navigators, and engineers closely identified with the practical aspects of 
iver work. Although obviously operating on all alluvial rivers, this phe- 
Omenon does not seem to be generally recognized among the hydraulicians, 


ote: Discussion open until January 1, 1958. Paper 1331 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 

_ Vol. 83, No. HY 4, August, 1957. 
Presented before American Society of Civil Engineers, at Annual Meeting, 

_ Jackson, Miss., on February 20, 1957. 

. Asst. to Chief af Eng. Div., New Orleans Dist., Corps of Engrs., U. Ss. 

_ Dept. of the Army, New Orleans, La. 

. Structural Engr., Bedell & Nelson, Cons. Engrs., New Orleans, La. 


1331-1 


e 


1331-2 HY 4 August, 1957 


and its effects have received little, if any, attention in evolving the theoreti- 
cal background of formal river hydraulics. 

As indicated, crossing bars, rise and fall in a manner broadly related to 
changes in discharge: On the lower Mississippi, the magnitude, of these 
changes ranges from 25%, to even as much as 50%, of the change in water 
surface level. This phenomenon is so sensitive that even a small rise is im- 
mediately reflected by a rise of the crossing bar. When a period of high flow 
is succeeded by a fairly rapid fall, these built-up crossings cause a river to 
be a series of pools with steep crossing slopes: Under the influence of the 
resulting high velocities, the crossings rapidly scour, with a resultant lower- 
ing of the pools to “normal” low water. Such scour may be disorderly, how- 
ever, and may fail to produce a good navigable channel; in fact, the more 
rapid the fall, the more likely it is that considerable dredging will be re- 
quired at the crossings. 

Closer study of this rise and fall of crossings, may throw new light on the 
movement of bed load material: It may develop that, in effect, the bends are 
“scoured out,” at high stages, and that much of the “scoured” material is 
simply moved to the first crossing. During lower stages, this material is 
probably gradually removed from the crossing and deposited in the bend be- 
low. In addition to this source of “back fill,” bends also receive back-fill 
from the extensive bank caving, which follows a falling river, after the sup- 
porting hydrostatic pressure is removed from the banks undermined at high 
stages. Consideration of these river bed “gyrations,” prompts a remark 
about hydrographic surveys; namely, that new bed configurations will be 
found just as often as such surveys are made, and the only really “complete” 
survey, would be a continuous motion picture of the river bed. 

This discussion of the rise and fall of crossings, in sympathy with changes 
in discharge, is presented for the purpose of relating it to the sand wave phe- 
nomenon, as one of the two major influences, affecting resistance to flow. It 
is hoped that the discussion will focus attention on a major hiatus in bed load 
studies; that is, the special nature of bed load movement in bends, which is so 
unlike bed load movement over crossings (or in glass sided troughs in labora- 
tories), that it may be considered as virtually a completely different type of 
movement; that is, there may be an intermittent movement from bend to 
crossing (or vice versa depending upon stage) and a continuous “through 
freight” movement, going on in parallel. 

The conclusion cited, concerning sand waves and sand wave systems, was 
influenced by several empirical sources: namely, (1) long experience in 
hydrographic surveys, originally by piano-wire lead-line and circle-sounding 
equipment, and, since the 1930’s, by one of the early applications of sonic 
soundings; (2) routine sailing-line fathometer channel patrols for navigation, 
(3) evidence, from a steamboat pilot house, of uniformly spaced “reefs” on _ 
the convex side of bends, having (at this location) gradual upstream slopes, 
over which the water flows at high velocity, and short abrupt downstream 
slopes, forming comparatively deep and still pools, which sharply check the 
flow, and (4) experience in paddling a canoe, rowing a skiff, and swimming 
among such “reefs,” when all these means of locomotion gave personal con- — 
firmation of the flow conditions described. 

The final conclusion is that these “reefs” extend into the channel, as deep- 
ly submerged sand waves, where they continue to produce the same effects as 
in the shoal water near the sand-bar shore; that is, they generate turbulence. 
This paper deals principally with these “reefs” as constituting a “primary” 
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wave system; that is, the “secondary,” and “tertiary” wave systems (both in 
evidence), are not given much consideration. The term “wave” will be re- 
stricted to these primary waves, and the secondary and tertiary waves will 
be referred to as “ripples.” 


Preliminary Field Investigations 


Continued speculation on the sand waves, led to the decision to make a 
preliminary field investigation: If the results warranted it, more extensive 
investigations could then be made, and any such preliminary work would 
determine the nature and direction of the more extended study. 

It should be noted here, that such investigations would be entirely imprac- 
_ tical except by the use of the fathometer to obtain a continuous profile along 
the path of travel; that is, no amount of conventional hydrographic surveying 
could yield the continuous detailed data necessary to study these sand waves. 
For this reason, previous “sand wave” studies, on the Mississippi, by 
McMath (around 1880) and Ockerson (about the turn of the century) are of 
little interest. 

A preliminary field investigation was arranged late in 1952: Advantage 
| was taken of a regular Navigation Channel Patrol trip of the stern-wheel 
Steamer General John Newton to obtain a continuous fathometer profile of the 
“sailing-line” from New Orleans to Old River, a distance of some 200 miles. 
| The horizontal scale ranged from 30 to 40 inches to the mile, and the vertical 
| scale was one inch to 10 feet. The observations, on this trip, confirmed the 
: general occurrence of sand waves over the crossings and around the bends. 


Hypothesis 


) The following 4-point hypothesis was formulated on the basis of accumu- 
| lated empirical experience and the results of the first fathometer “sailing- 
line” profile. 

1 Point (1) 

| Sand waves occur generally on the bed of a river in a systematic manner, 
; and they vary in a systematic manner with changes in discharge. 

Point (2) 

| The sand waves tend to become larger with increases in discharge and 

i smaller with decreases in discharge. 

| Point (3) 

There is a lag in the rate of change of sand waves, and, on a continuously 


rising river, the waves are somewhat “smaller than they should be,” and, on 
a falling river, they are somewhat “larger than they should be.” 


Point (4) 
These waves constitute a major element of resistance to flow. 


Corollaries 
_ These corollaries proceed from the hypothesis: 


Corollary (1) 

Roughness, or resistance to flow, is not an accident or an attribute of the 
soil formations through which a river flows, but is rather a property, or 
-ondition, of the river bed, controlled by the river itself. 
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Corollary (2) ! 

The fact that the “rating curve,” at any given station, is a “loop” (Fig. 1) 
rather than a “curve” may be partially explained on the basis of a “lag” 
sand wave adjustment; that is, the smaller waves present, on a rising phase, 
offer less resistance to flow, and the larger waves present, on a falling phase, 
offer more resistance to flow. 

The over-lapping phenomenon, of rise and fall of river crossings, also 
operates to produce the conditions recorded by the “loop” rating curve. 


Corollary (3) 

Slope itself (in a local sense) tends to become a “result” rather than a 
“cause.” 

Underlying all this discussion, is the assumption that, with discharge and 
other governing conditions constant, a stream will produce a system of sand 
waves, which remains constant, as to size and shape of individual waves, and 
as to the wave “system.” The implication is that, during a cycle of rise and 
fall, bed adjustment cannot take place, with sufficient rapidity, to be, at all 
times, in accord with the instantaneous demands of the varying discharge. 


Extended Field Observations 


Since the initial field observation in 1952, observations, on an extended 
basis, have fully confirmed that sand waves exist in a systematic manner, and 
vary in accordance with some orderly influence. Although the “systematic” 
occurrence is clear from visual inspection of the profiles, the influences 
governing the “system” are less clear. There are, however, definite indica- 
tions that “individual waves” and “wave systems” vary in sympathy with 
changes in discharge, and there is some likelihood that a certain type of wave 
and wave system will be found to be associated with the horizontal alignment 
at any location. 

As to the effect of increased discharge in increasing wave size, there was 
early evidence of an “anomaly;” that is, it was noted on two occasions, when 
the river rose to nearly bankful, that the “tops” of some crossings became 
perfectly smooth: The word “anomaly” is employed because waves were 
still present, on the upstream and downstream slopes of the crossings and 
around the bends, and the waves, where present, were larger than at low 
water. This behavior, of the tops of crossings, will be an interesting facet 
of study to pursue further, when high over-bank river stages again occur. 

As a first result of these field observations, the following tentative deduc- — 
tions are advanced: The direct effect on resistance to flow, of the size and 
shape of the particles on the river bed (i.e., those on the immediate surface 
at a given instant), may be small in comparison with the effects of the suc- 
cessive orders of waves and ripples. Since the “primary” sand waves are so 
large, as to virtually constitute major river-bed “topography,” it appears 
possible that their frictional effect is great as compared with, either the ef- 
fect contributed by the lesser orders of waves (i.e., “ripples”), or the 3 
theoretical “drag” effect. Finally, it may be that the sizes and shapes, of . 
grains, are ultimately significant, not because of their direct effect upon fric- 
tion, but because of their indirect effect, in governing the formation, size and 
systematic variation of the primary sand waves and “ripples” on the river 
bed. It also appears likely that the size range of particles, present onthe — 
river-bed surface, at any given time, is subject to control by the river, oe 
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as seems to be the trend of opinion in recent papers, based on laboratory and 
field measurements, and analytical studies. 

The observations, on which this paper is based, were entirely fathometer 
profiles; that is, no studies were made of sand sizes, the manner of adjust- 
ment of the waves, their speed of movement, or other associated phenomena. 


Views on Sources of Friction in Alluvial Rivers 


Various hydraulic texts and published papers, treating flow of rivers, bed 
load movement, and related questions, have been examined in order to 
develop the thinking with respect to systematic variations in the river bed 
(i.e., sand waves, and the rise and fall of crossings) as major phenomena 
governing resistance to flow and the general behavior of alluvial rivers. In 
general, the broad subject, of such systematic changes, does not seem to en- 
ter into formal river hydraulics: Rather, there is extended consideration of 
friction as such, and a change in attitude, over the years, as to the causative 
factors. 

The generally accepted view, of earlier writers, appears to have been that 
“roughness” is a property, inherent in the channel, and ultimately related to 
the nature of the soils, or the geological formations, through which the river 
carves its channel. The impression seems to have been that, whatever it was 
that controlled friction, it was constant for a given reach, or even an entire 
river, once its value had been determined. 

Another group of writers, relates the friction to the “sand paper” or “drag” 
effect of the particles, present on the surface of the bed. Some go further and 
indicate a belief that the river exercises a selective influence over the sizes 
of grains, actually present under any particular condition. Still other writers 
take a position, somewhat related to that presented in this paper; that is, sand 
waves and ripples, are a factor in determining resistance to flow. At least 
one has indicated the possibility that the waves get larger with increased dis- 
charge. Another thought the river bed would be smooth at the lowest stage. 

Still other writers make river bed behavior very complex: They put un- 
limited trust in involved mathematics which finally rests on field observa- 
tions, or laboratory measurements, that may contain an element of dubious- 
ness. By insisting on a detailed particularized approach, they implicitly 
reject the possibility of a generalized approach, based on laws or principles, 
perhaps susceptible of initial identification and subsequent measurement on — 
an empirical basis. It may be, figuaratively speaking, that microscopic work 
on sand grains, and small areas of stream beds, has prevented recognition 
of (1) the need to view the river bed continuously in profile; (2) the possible 
significance of the systematic river bed changes discussed in the opening | 
paragraph of this paper, and (3) the possibility that bed-load movement in the - 
bends, is so different from bed-load movement in the crossings, that they 
constitute virtually separate and distinct phenomena. : 

Summing up, it appears that the proposed hypothesis conflicts somewhat i 
with both the older and the more recent prevailing views concerning causes 4 
of flow resistance, by implying that the major source of resistance arises 
from systematic river-bed changes, which the river makes, of its own voli- 
tion, and modifies to meet changing discharge requirements. The overall 
conception is a river which lets itself down its valley under control, much as 4 
a person controls his descent, down a flight of stairs, with some additional 
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control obtained by sliding his hand down the bannister. Final support, or 
| partial refutation, of the four-point hypothesis, awaits additional fathometer 
| observations at overbank stages, and comparisons of the resulting profiles 
with similar profiles, taken at lower stages. Because of the prolonged low 
stages, due to the drouth, no overbank stages have occurred in the five years, 
Since interest in this subject came to a head. 


Comparisons of Fathometer Sailing-Line Profiles 


Successive Sailing-line profiles were run, over several years, under vari- 

| ous conditions of discharge. These profiles have been studied to detect sys- 
tematic variations, and to select half-mile “study-reaches” for comparing 

such variations. The results of these studies have been incorporated in the 

figures, now to be discussed: 

| Each figure has a vicinity map to show the river configuration above and 

below the selected “study-reach.” There is also a condensed “sailing-line” 

| profile, extending above and below the study-reach. The location of the study- 

_ reach itself is shown on the map and on the profile. Finally, the half-mile 
study reach itself is shown to full scale, for each of the conditions, to be 

| compared. 

In viewing the full scale study-reach profiles, it should be noted that the 
Sand waves are somewhat distorted: This is due to the fact that the “vertical” 
axis of the fathometer chart is a curve; hence, it is “vertical” only at the cen- 
ter of the chart, is inclined up and to the left on the lower half, and is inclined 
up and to the right in the upper half of the chart. The vertical scale is as in- 
dicated by the depth numerals, and the horizontal scale is as indicated by the 
marks, which set off 20ths of a mile, or roughly 250 feet. 


| Figs. 2, 3, and 4 

A fairly typical crossing: In fig. 2, at the higher stage, the upstream slope 
. has smoothed off, but there is a well developed wave system on the rest of the 
| Crossing, and the waves are larger at the high stage than at the low stage. 
Note that the amplitude of the waves, on figs. 3 and 4, is nearly 30 feet and 

' that the normal primary waves have formed on these “Super” waves. It is 

| Suspected that the “Super” waves, just forming here, will take over complete - 
| ly at higher than bankful stages. In fig. 4, the waves, on the downstream 

/ crossing slope, are of about the same amplitude at both stages. This devia- 

| tion from the hypothesis, cited earlier, does not affect its general validity, as 
| will appear in due course. 


| Figs. 5 and 6 

Both of these figures deal with the lower end of a crossing at the start of a 
bend. In fig. 5, the waves on the higher stage are larger than the waves on 
the lower stage. In fig. 6, the waves are much larger at the higher stage than 


| at the lower stage. 


Figs. 7, 8, 9, and 10 
To the right of fig. 7, the Sand Wave System for the higher stage, has built 
up to an amplitude of about 10 feet. These waves are not present at the lower 
stage. Note that the crossing has smoothed off at the higher stage. In fig. 8, 
a “Super” Sand Wave System has attained an amplitude of over 20 feet, at the 
higher stage. Note that the normal primary waves are superimposed on the 
“Super” system. At the lower stage there is no significant wave system. In 


1a 


HY 4 August, 1957 


8 


1331- 


O'l9| OL G’9si SATIN 


YSAIN IddISSISSIW 


SNIT-ONITIVS = 


a 


-9 


1331 


KELLER 


CAREY 


CE 


AS 


a aa 


ok 


ee ee 


9S6l Wddv Ol 3WsOdd J3NIT- 


eos l 


YaAIY 


4@S ose 


Sa tier cece e c 


pn 


OL O 


An 


191 


Fie 


poor. 


* 


Ee 


et 


SW3LSAS 3AVM GNVS _ 


August, 1957 


4 


1331-10 


09 mo . Nemriiee taccr es er S 
ls ae 8! Q9S61 Tiddv O| a6e9 “BIIAUOSPIOUOG 


ns  SWALSAS SAVM GNVS— 


SM mojag 1287 syceg 


9S6l lddV Ol AlNsOdd ANITI-SNIMIVS 


€9S| OL O'19I ATIW 


YAAIa ime 


SNIT-ONIMIVS — 


<= a | Liaise res 


HOVSa AGCNLS 


1331-11 


CAREY - KELLER 


r 


— 9S6l acy ol 


ee : 


LEASE GASES 


G°cer 


O'Sel OL O'GEL ATIIW 


daAle iddiSSissim 


, 19 


August 


4 


HY 


-12 


1331 


Bo lee Mate Ye oe oe 
™ shod SIASR OT cy 


Atel og6rtddv Ob 


9S6l 1lddv Ol AWAOYd ANIT- ONIiVvs 
O'SZ} OL O'SEl 3TIW 


YsSAla Idd ISSISSIN 


1331-13 


CAREY - KELLER 


sso 


 NNgpted@S S61 AINE Ob sages 


See 


9p aijiAuospjouog Z°991 
oot sie i gen A3€e°Si 9S6I Wadv Ol E 9g} 


SWALSAS ZAVM 


—_ WHEL vi czeczzaa. 


S°vel Sob S°seL 98h s°9el 48h 


9S6l 1lddv Ol AN4ONd ANIT-ONMIVS | | 


B64 OL S48 STIW 


YAAIM IddISSISSIN 


eae 
\ } 3NIT-ONITIVWS 


SH mojieg 4284 ‘sydd 


August, 1957 


HY 4 


1331-14 


rvol i abp9 9||[AUOSpjou0g : : : Sah : 


oo} Bae nies Pe 


991 s°9el 48} 


9S6l 1lddv Ol JWsOdd ANIT-SNIMIVS 


S"°6Zl OL S2Z8l JTIW 


ere lacs S12 a 


SM Molag 1984 :syscog 


1331-15 


CAREY - KELLER 


Op aa ee 
- . 39S 9s 9s6l mar OL ae 


abp9, a}|!Auospjpu0g eel 


aa SERIE ieee 


BEES ERM: Ol. 


SWNALSAS JAAVM GNVvS te 


oO 


S‘1el Zel 


9S6I WddVOl 31ldOudd ANIT- ONIMIVS 


8’6Lh OL SB SATIW 


feet tte Jd¢|2 oc 


195 


? 


August 


4 


HY 


1331 


Ovl 
mete: 
O92 


= 
> OSI gt S*Iel eel 


9S6I TIddV O1 


S*zel 


31140ud ANIT=9NITN 
S6lt OL S%Sl SJTIIW 


haere el Id ISSISSIN , 


ANIT-ONITIVWS 


zai! ae 


HOVad Ads 


ASCE CAREY - KELLER 1331-17 


ig. 9, the explanation for the system of large waves, at the lower stage, is 

t this is the same system of large waves, as shown in fig. 8, moved down- 
tream. As the stage drops, it is possible that the waves migrate down- 
tream, decreasing in size as they move. In fig. 10, the waves in these two 
tystems are not too different; possibly because of the equalizing effect of the 
ave migration, as indicated in figs. 8 and 9. Note that the downstream 
rossing slope has been somewhat smoothed off, at the higher stage. 


ig. 11 

‘This figure shows “Grandview Reach,” a long “reach” that has been his- 
orically stable. At the far left, for the highest stage, there is a well 
eveloped wave system with an amplitude of about 10 feet; at the middle 
tage, the wave system has migrated to the center of the study reach and the 
mplitude has decreased to about 4 feet; and, at the lowest stage, the major 
= system has migrated further downstream and the amplitude has de- 


reased to 2 or 3 feet. 


ig. 12 

This figure shows a rapidly caving bend above Baton Rouge: This bend is 
omposed almost entirely of sand to the top of the bank. The unusual shape 

the waves at the high stage is due to the reduced horizontal scale of this 
igure. The steep slope from the lower end of the crossing to the head of the 
end has smoothed off at the high stage. This bend has caved at the average 
ate of 5,000,000 cubic yards per year for the last 16 years and has moved 
ack 70 feet per year at the point of greatest caving. 


igs. 13 and 14 
Note the unusual alignment of the river in this vicinity. The waves are 
ger on the higher stage. In fig. 14, there are indications of “Super” waves, 
which the normal primary waves are superimposed. 


igs. 15 and 16 

_ These figures show a typical crossing: Fig. 15 shows the upstream cross- 
slope, and fig. 16 the downstream crossing slope. In fig. 15, there is a 

ell developed “Super” wave system, on which are superimposed the normal 
rimary waves. However, at the lower stage, only the normal primary waves 
re present. Fig. 16, shows rather symmetrical and consistent wave systems, 
t both stages. There are also indications that there has been an “equaliza- 
.on” of the wave systems, at higher and lower stages, due to the migration of 
ne larger waves from the upstream crossing slope. 


| CONCLUSIONS 

These studies afford substantial support for three of the component points 

ee hypothesis; namely, (1) that sand waves exist generally on the river bed 

a systematic manner; (2) that sand waves tend to vary in amplitude and 
acing, in accordance with variations in discharge; and (3) that a major 

is of resistance to flow, in an alluvial river, is provided by a primary 

ystem of sand waves, which, in turn, is created by the river itself and oper- 

‘eS as a means of controlling its orderly progression down its valley. The 

urth “point,” that the waves constitute a major component of resistance to 

ow, has not been confirmed, although it still appears reasonable. 

: The fathometer profile observations, on which this paper is based, covered 
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a limited range of river stages, and the value of the conclusions, derived 
therefrom, is correspondingly limited. The paper will serve its purpose if i$ 
demonstrates the existence of a practical field of investigation that may lead 
to useful knowledge concerning; (1) the adjustment of sand wave formations, 
(2) the rise and fall of crossings, (3) the special character of bed load move- 
ment in bends, and (4) the joint bearing which these phenomena have upon re- 
sistance to flow. 

Federal Agencies, engaging in planning, constructing and operating water 
resources projects, can often gather data on river bed conditions at little ex- 
pense over that of the fathometer paper involved; as, for instance, in the case 
of channel patrols for navigation. In general, observations can be made at 
full cruising speeds on trips for other purposes, or on the idle return trips 
from other missions. Either the “sailing-line,” the “thalweg,” or the river 
“center-line,” can be followed, as may be convenient, and it is not necessary 
to cover an “area” in order for the data to be useful. Working up the data, 
can conceivably be relatively inexpensive. A point to note is that studying 
“lines” or “profiles,” has much to be said for it, as against studying “cross- 
sections”: In the first place, the key to studying these river bed phenomena 
is “continuity,” and this is something “cross-sections” cannot supply, despite 
the great amount of field and office work involved; in the second place, a 
“line” can be taken as representative of a “volume” within certain limits, if 
proper lateral-control of the line is provided, and known points are selected 
for linear-control “fixes.” 

After a certain amount of “qualitative” study of “lines,” to determine the 
manner and magnitude of river bed adjustments, ways and means will prob- 
ably be evolved to use the same data to make “quantitative” determination of 
the bed-load material movement involved in the sand waves and in the rise 
and fall of crossings; that is, the material that probably has the most pro- ~ 
found effect on the hydraulic properties of streams. 

In conclusion, it appears likely that organized study of these phenomena of 
river bed behavior, might have these effects: (1) throw new light on some 
theoretical aspects of river hydraulics, (2) indicate a better approach to plan 
ning and constructing bridge piers, pipe-line crossings, etc., and (3) point the 
way to a better understanding of the hydraulic results of comprehensive rive 
stabilization and lay the basis for extending the claimed economic benefits, 
for such work, far beyond the mere protection of lands and improvements 
from the effects of bank caving and bar building. 
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SYNTHETIC STORM PATTERN FOR DRAINAGE DESIGN 


Clint J. Keiferl and Henry Hsien Chu, 2 Associate Members ASCE 
(Proc. Paper 1332) 


| SYNOPSIS 


The City of Chicago has been engaged in an extensive auxiliary outlet 

| Sewer construction program since the end of World War I. Bond issues 
totaling more than $88,000,000.00 have been passed to finance this program. 
| The engineers of the Sewer Planning Division felt it worthwhile to make a 

| comprehensive study of the principles involved in the rainfall-runoff relation- 
_ Ship in order to economically utilize this expenditure. The development of a 

_ Synthetic storm pattern, for use in the hydrograph method of sewer design, is 
one phase of this research, and is the theme of this paper. 


INTRODUCTION 


It is the purpose of this paper to present a method of determining a storm 
pattern which may be used in the hydrograph method of sewer design. (1,2,3) 
This storm pattern will be developed from a rate-duration-frequency curve 
Selected for use in design. It will also be shown that although this storm pat- 
tern will include the average intensities of the rate-duration curve for all 
durations, it will produce no greater runoff peaks than that of the rainfalls 
Of the separate durations. 
| The chronological location of the peak period of rainfalls with reference to 
the total storms period and the amount of antecedent precipitation immediate- 
ly preceding the maximum period of any durations will be derived from the 
Statistical average of rainfall records. The importance of this antecedent 
rainfall in affecting the peak of the runoff hydrographs will be illustrated. 


Note: Discussion open until January 1, 1958. Paper 1332 is part of the copyrighted 
_ Journal of the Hydraulics Division of the American Society of Civil Engineers, 

_ Vol. 83, No. HY 4, August, 1957. 

1. Senior Sewer Designing Engr., City Bureau of Eng., Chicago, Il. 

2. Civ. Engr., City Dept. of Public Works, Chicago, Il. 
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Rainfall Characteristics Affecting Runoff 


The rainfall and resulting runoff hydrograph should be analyzed in order 
to determine what other characteristics, besides the average rate of water 
falling within a given time period, would tend to affect the peak rate of run- 
off. If in one rainfall, some precipitation preceeded the maximum period, 
filling many of the surface depressions, wetting all surfaces and satisfying 
the greater infiltration capacity of the soil, the peak runoff rate would be 
higher than from another rainfall having the same intensity during the maxi- 
mum period but lacking the antecedent rainfall. In the development of a syn- 
thetic storm pattern for use in the hydrograph method, the right amount of 
antecedent rainfall must be replaced in front of the statistical volume of rain- 
fall occurring within the maximum period. 

Another important characteristic that affects the peak rate of runoff with a 
given volume of water falling in a specific period, is that of the variation of 
rainfall intensities within that time period. Consider a long duration rainfall 
in which most of the precipitation falls in the early part of the period. In this 
case, the greater infiltration capacity at the beginning and the surface depres- 
sion storage would absorb much of the peak rainfall, thus making the runoff 
rate approach a lower value of equalibrium. On the other hand, if the bulk of 
precipitation occurred in the latter portion of the period, most of the previ- 
ously mentioned losses would be already satisfied before the time of peak 
rainfall intensity, and a higher peak rate of runoff would result. 

To sum up briefly, the three most important characteristics affecting the 
peak runoff rate for a specific period or duration are as follows: 


1. Volume of water falling within the maximum period, 
2. Amount of antecedent rainfall, 
3. Location of the peak rainfall intensity. 


A synthetic rainfall encompassing the statistical average of these three char- 
acteristics should then be adequate for use in the hydrograph method as the 
design storm pattern for any given duration, see Figure (1). 


Evaluation of the First Characteristic for Single Storm Pattern 


The volume of water falling within the maximum period can be taken from 
the rate-duration curve of a given frequency selected for design purposes. 
The equation for this rate-duration curve may take the following form: 


RM ali = Sr) 


where igy, is the average intensity in inches per hour, tg is the duration of the 
maximum period in minutes, and a, b and c are constants. The volume of 
water falling, or mass rainfall for any duration is graphically represented by 
the area of a rectangle of which the upper right-hand corner touches the rate- 
duration curve, see Figure (2). This rainfall mass, in inches, can then be ex- 
pressed by 
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Substituting Equation (1) for igy 


pet eee 
eid Fomeary) (2) 


First, a completely advanced type rainfall pattern can be developed which has 
the peak occurring at the beginning, and with no antecedent rainfall. A curve 
representing the hyetograph of this type of rainfall must have the same aver- 
age intensity for any maximum period as that given by the rate-duration 
curve, as shown in Figure (2). The area under the curve of the hyetograph 
can be expressed mathematically by 


4 
P= abl vat 
0 


where i is the ordinate of the hyetograph in inches per hour. 
Differentiating with respect to tq 


| df Pie a 3h (3) 
PTA, eters 
| But from Equation (2) 
Sate Se bl? 
fae & CO ty eee 


| Differentiating 


| dp _ g [l-bi+c] (4) 
| dt, 60 btea.c)” 


| Combining Equations (3) and (4) 


b 
ris a (i-b)t i c] (5) 


| x GS + c)e 


| Equation (5) then represents a completely advanced type storm having for all 

durations, the same average intensities as those given by the rate-duration 
curve. For Chicago area, the accepted rate-duration curve developed by 
Eltinge and Towne of the Sanitary District of Chicago, (4) for a five-year fre- 


quency, is 


90 (6) 
Chay £994 WW 


Therefore, the constants in Equation (1) are a = 90, b = 0.9, and c = 11. 
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Substituting these values into Equation (5), one obtains 


90 (01 6°74 11) (7) 
(477+ I1)* 


By assuming different values for tq, Equation (7) represents a storm starting 
with the rate of 8.18 inches per hour at the first instant and decreasing 
steeply to 4.42 inches per hour at the end of the 5th minute, to 1.13 inches 
per hour at the end of the first half hour and so on. 

To make Equation (5) applicable to any intermediate type storm pattern 
the following modifications shall be made: Within the maximum period of 
any rainfall, the duration tg can be split up into that part occurring before 
the most intense moment and that part after the most intense moment. In- 
troduce the symbol ‘r” to represent that portion of any duration occurring 
before the most intense moment, expressed as a ratio to the entire duration, 
refer to Figure (3), and set 


Z=rk and = d-nt (8) 


where tp is the time before the peak in minutes, and measured from the peak 
to the left and tg is the time after the peak in minutes measured to the right 
of the peak. The ratio “r” is then a measure of the advanceness of the storm 
pattern. If r = 0, the storm is of the completely advanced type as discussed 
in the last paragraph. If r = 1, the storm is of the completely delayed type, 
that is, it has its peak at the end of each and every duration and has con- 
siderable antecedent rain falling before each and every maximum period. 
For intermediate type storm patterns, the “r” value will range between zero 
and one. 

Solving the first part of Equation (8) for tq and substituting in Equation (5) 


rye 
(Before The Peak) fn; ve a{(-) (4) z c] (9) 
[(B)° + <]? 
sea the second part of Equation (8) for tq and substituting in Equation 
5 
a[(-(5)"c] | 
(al J a 


fav 


(After The Peak) 


A synthetic hyetograph plotted from Equations (9) and (10) will have, for 
all durations taken during the most intense period, the same average intensity 
as the rate-duration curve from which the constants a, b and c are derived. — 
Using the above mentioned constants for the Chicago area, 
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90 (0.1 (72 yu) 
(After The Peak) hee — a . (12) 
[()°% 11) 
Determination of “r” Conforming with the Second and Third Characteristics 


The first characteristic has been evaluated by the development of Equa- 
tions (11) and (12) from the accepted rate-duration curve. It remains to de- 
termine, from records of excessive rainfalls, the best value to assign the 
constant “r” so that the other two characteristics of the synthetic storm pat- 
tern will conform approximately to the statistical data. 

Rainfall records were obtained from the Metropolitan Sanitary District of 
Greater Chicago, where the records of heavy storms from seventeen automa- 
tic recording rainfall gauges located in the Chicago area are compiled. 

Four widely separated stations were selected so as to retain the individuality 
of each station-rainfall. The four stations were at Loyola University (Station 
10), Springfield Avenue Pumping Station (Station 14), The West Side Treatment 
Works (Station 7) and Roseland Pumping Station (Station 15). Only those 
rainfalls which at least for some duration are considered excessive, were re- 
corded. Eighty-three station-rainfalls are given, with data taken from the 
records for durations of 15, 30, 60, and 120 minutes. 

In Table 1, the following information is given: For example, at Station 14 
on May 11-12, 1935, the maximum rainfall occurring in a fifteen minute 
period was 0.39 of an inch, while the antecedent rainfall occurring before that 
period was 1.40 inches. Since the maximum fifteen minute period contains 
three five-minute periods, the most intense was the second. Likewise, in the 
same storm at Station 14, the maximum mass occurring in the 120 minute 
duration was 1.82 inches, only 0.05 inches fell antecedent, and out of the 
120 minute duration the twentieth five-minute period was the most intense, 
therefore the storm was of the delayed type pattern. The blanked spaces indi- 
cate those durations within which the maximum mass rainfall was not exces- 
sive. These figures are purposely deleted in order that the statistical aver- 
age obtained will reflect the characteristics of only the excessive storms 
concerned. The mean values of antecedent rainfall and the location of the 
peak for each of the given durations were obtained. The mean values of mass 
antecedent rainfall are plotted in Figure (4), and the mean values for the loca- 
tions of peak are plotted in Figure (5). | 

An analytical expression for the antecedent rainfall of the synthetic storm — 
pattern can be derived as follows: For a given “r” value, the mass rainfall 
within the duration tq occurring before the peak can be obtained by multiply- 
ing Equation (2) by “r,” or 


where Pp is the mass rainfall before the peak in inches. _ ¥ 
For the Chicago five-year frequency rainfall 


Pp, = £L5 l= 
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Then the total duration of the synthetic storm pattern is set equal to the long- 
est time of concentration of any sewer system to be designed. For Chicago 
the total duration of the synthetic storm pattern was set at 180 minutes. Also, 
since this duration was quite long for storm sewer design, it was assumed 
that any small amount of rainfall preceding this period would produce little 
fis 1 on the rate of runoff, therefore no antecedent rain was assumed before 

Ss 180-minutes-long synthetic storm pattern; then the amount of antecedent 
rainfall preceding the duration tg, can be expressed as 


: 
: A =[R] CR] = r ee rt5@ ~rh29- ete mee) 


“= 180 £94 I/ ZG 9 


here A is the mass antecedent rainfall in inches. Solving Equation (13) for 
r” and substituting the statistical mean values of antecedent rainfall for A 
dd the corresponding duration tg, the value of “r” for each duration can be 
computed, see Table IL. 


| 1.000 |1.290 
| {= 


1.588 | .902 


| .8 | 1.769 | .521 | .203 


4 


| 


Bello | Lso 118 


Table II - Determination of the Value of“r” from Mass 
Antecedent Rainfall. 


lince the “r” is to be constant through all durations, and the antecedent rain- 
ll for the shorter durations is much greater and of greater importance in 
itfecting runoff than that for the longer durations, it seems that the “r” com- 

ited for different durations should be weighted in proportion to the mass it 
eflects in obtaining the average “r” value. The weighted average value of 
ip thus obtained is 0.386. 

The location of peak within each duration provides another direct approach 

\ evaluating the “r” value. In the last line of Table I, “Mean Location of Peak 
Rainfall Duration” or mean “r” for each duration, is calculated on the 
is that the momentary peak of every recorded rainfall happened exactly at 
e middle of the recorded “Peak Five-Minute.” Then these “r” values of 
ch duration were weighted proportionally to the time of duration in obtaining 
e average “r” for all durations. This was done because the shorter dura- 
ms contain so few five-minute periods that the mean location of the peak is 
So well defined as in the cases of longer durations. The value of “r” thus 
rived by considering the location of the peak is 
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[5 (453) + 30(4/5) + _60 (356) + 1206369 _ 9 376 
/§ + 30 + 60+/20 


An even fraction 3/8 or 0.375 which is close enough to the above findings 
was assigned to “r,” thus establishing a synthetic storm pattern with a three- 
eighths advanced type. The mass antecedent rainfall for “r” = 0.375 is plot- 
ted in Figure (4) by the use of Equation (13). Also in Figure (5) a horizontal 
line with “r” = 0.375 is drawn. Comparison of these curves with the statisti- 
cal mean values indicate fairly close agreement in antecedent rainfall with 
the shorter durations and in the location of peak with the longer durations. 

Substituting the value of “r” = 3/8 into Equations (11) and (12) they become 


8 0.9 
(Before the peak) ; es 90 te * + # 
SC) etl 


ao[oi(ot,)”* i] 
(84)? f° 


The synthetic storm pattern plotted according to the above equations satis: 
fies best the three characteristics outlined in the beginning of this treatise 
and therefore can be used as the design storm pattern for the hydrograph 
method of storm sewer design in the Chicago area. Since the total duration 
of the synthetic storm pattern is limited to 180 minutes, the time before the — 
peak is limited to 3/8 of 180 minutes or 67 1/2 minutes and the time after the 
peak to 5/8 of 180 minutes or 112 1/2 minutes. A hyetograph of the synthe 
storm pattern is shown in Figure (6). 

Also shown for comparison, in Figure (6), is a hyetograph for a one-year 
frequency synthetic rainfall, based upon the Eltinge-Towne formula for a 
“One-year frequency” rate-duration curve, namely 


(After the peak) fetes 


ee = ae 
eV. Pipe 7 


and the value of “r” is assumed to remain at 0.375. 


Pe eae > Re he = 


CONCLUSION 


The “Synthetic Storm Pattern,” having been derived from the accepted 
rate-duration curve, and fitted to the statistical values of antecedent rainf. 
and location of the peak, may now be used in the “Hydrograph Method” of 
sewer design. The advanced type pattern with “r” = 3/8, found to be accept: 
able in the Chicago area, might not be suitable in other localities. The au- 
thors believe that setting the “r” value is of considerable importance since 
it determines the amount of antecedent rainfall for all durations. In Figure 
(7), it is shown for comparison, the affect of the change in the “r” value on 
the runoff hydrograph of a 70 foot length of turf strip after subtracting losse 
due to infiltration, depression storage and surface detention. It can be seen 
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at for the completely advanced storm pattern, the peak runoff rate is only a 
ird of that for the 3/8 advanced storm pattern. Figure (8) shows a compari- 
F of the runoff hydrographs from a typical residential area in Chicago, 

th 36 percent imperviousness, resulting from synthetic storm patterns hav- 
' er” = 0.375 and “r* = 0. 

The Synthetic Storm Pattern having the same average intensities as that 
ven by the rate-duration curve for all durations, would not be likely to occur 
the same frequency as the rate-duration curve from which it was derived. 
would undoubtedly have a greater return period, that is, less frequent occur- 
nce. The use of the Synthetic Storm Pattern, however, should not produce 
eater runoff peaks then the use of a separate uniform intensity rainfall for 
ery location in the sewer system, if the uniform intensity rainfalls have the 
rrect antecedent precipitation. In Figure (9), the synthetic storm pattern 
S been used as a supply hydrograph at intervals along a sewer, and routed 

the time-offset procedure to produce hydrographs at various times of 
avel. This use of the synthetic storm pattern, and subsequent use of uni- 

m rainfalls as the supply hydrographs, is only illustrative and would occur 

y if the area was entirely impervious and had no losses or reductions be- 

e entering the sewer. Figures (10) and (11) show various duration uniform 
‘ensity rainfalls, with antecedent precipitation equal to that of the 3/8 ad- 

ced synthetic storm pattern, being routed by time-offset to produce hydro- 
.* at various times of travel. Figure (12) shows the comparison of the 
s of the hydrograph derived from the synthetic storm pattern and those 
ived from the uniform intensity rainfalls. It can be seen that the attenua- 
in of the hydrograph peaks of the synthetic storm pattern envelope those of 
> uniform rainfalls of separate durations, and would produce no greater 
tities for use in design. In a more rigorous analysis of comparing the 
ign rates derived from the synthetic storm pattern to those derived from 
infalls of separate durations, it was found that the peak rate of runoff is 
iy proportional to the net volume of the supply curve in the maximum 

iod equal to the routing time. Then, since the synthetic storm pattern 

the same volume of rainfall up to the end of any duration as the uniform 

msity rainfall of the same duration with antecedent precipitation, they 

Jd both produce the same runoff peaks. If infiltration, depression storage 

‘surface detention were subtracted chronologically from these two rain- 

s, the same net volume of supply would result. Thus after routing, they 
uld again produce equal runoff rates. 
iFrom the foregone, it is concluded that the synthetic storm pattern can be 
ed as the “Design Storm Pattern” in the hydrograph method of sewer design 
, all locations in a sewer system. 
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SYNTHETIC STORM PATTERN EMCOMPASSING THE 
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_ CHARACTERISTICS OF FLOW OVER TERMINAL WEIRS AND SILLS 
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| ABSTRACT 


| Generalized experimental results are presented to show the variation in 
ith the discharge coefficient and the nappe profile for two-dimensional flow 
yer a vertical sharp-crested weir at the end of a horizontal channel as the 
itio of head to depth of flow changes continuously from zero to unity. 


| SYNOPSIS 


| Whereas most existing weir formulas lose significance as the ratio of head 
jheight of weir becomes great, it is known that the rate of flow over a sill at 
2 end of a channel can be estimated from the basic weir equation by assum- 
the velocity of approach to have its critical magnitude. This gives rise to 
ischarge function for sills which meets that for weirs in a sharp peak when 
» head is about ten times the weir height. There are related indications 

t the nappe attains a maximum elevation and thickness at the same head- 
ight ratio. Laboratory tests to determine quantitatively the variation in 
scharge coefficient and nappe profile over the entire range are described, 
a the experimental results are presented in generalized form for design 


poses. 


Historical Background 


Poleni,(1) in 1717, was the first to publish an equation for weir discharge 
ained by integrating over the outlet section an expression for the velocity 
rerms of the square root of distance below the free surface. Although he 


: Discussion open until January 1, 1958. Paper 1345 is part of the copyrighted 
ournal of the Hydraulic Engineering Division of the American Society of Civil 
mgineers, Vol. 83, No. HY 4, August, 1957. 
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neither took account of the nappe contraction nor explicitly introduced the 
gravitational acceleration g, Poleni’s name continues to be associated with 
the elementary discharge equation 


q= S C Vee 23/2 (1) 


in which q is discharge per unit width of flow section, C¢ the coefficient of 
contraction, and h the head on the crest. Contraction of the nappe was taken 
indirectly into consideration by Du Buat(2) sixty years later in setting as ar- 
bitrary limits of integration the crest level and a level the distance h/2 above 
the crest; it was also Du Buat who introduced present-day methods of estimat 
ing the effect of submergence. Bidone, (3) better known for his study of the 
hydraulic jump, appears to have been the first to show an interest in express- 
ing algebraically the geometry of the weir nappe. 

Not until about the middle of the 19th century was attention focused upon 
the velocity-of-approach “correction” for weir discharge, which was to domi- 
nate certain branches of the literature for another century. In 1845 
Weisbach(4) published the familiar integral form 


3/2 2\3/2 
2 v (2) 
Fev [b+2) -(2) | 


Since the solution of this equation has to proceed by successive approximatior 
because of the implicit presence of the approach velocity vg on the left side 

of the equation, the impression has persisted that this velocity is independent 
ly variable. Weisbach himself appreciated the fact that a unique solution . 
existed for any ratio of head h to depth of approach h + w, as he proposed an — 
empirical formula for its evaluation: 


q= 


(c.}5 [1-08 + 0.3693 (sis) [var (3) 


wiloo 


Bazin(5) in 1888 adopted the same form of function, with specific evaluation 0 
the limiting coefficient of contraction (Cag 4 


a= (0.405 + 9:9°5)[ 1 4 0.55 (2 st) | va 0? 4 


the dimensional term 0.003/h being introduced to compensate for low-head € 
fects now attributed to surface tension and viscosity. In a series of investi 
tions extending over the next ten years, Bazin also went deeply into such 
varied aspects of weir flow as the distribution of velocity and pressure 
through the nappe, and the effects of incomplete ventilation, submergence, ~ 
weir inclination, and geometry of the crest; his nappe-profile measures 
are still used as the basis of spillway design. d 
Bazin’s formulation of the discharge coefficient soon led others to atte 0 
further improvements by slight changes in the type of algebraic expression 
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ised. From 1911 to 1928 Rehbock(6) proposed a number of variants, the 
nest pertinent of which was 


h 3/2 
q= (0.60s + 0.08 aes ERE EE gh / (5) 

: 
fomparison of Eq. (5) with Eq. (4) will show that the primary distinction lies 

the use of the ratio h/w in place of h/(h+w) as the geometric parameter. 
he remaining term of Eq. (5)—which, like the corresponding one of Bazin, 

not dimensionally homogeneous with the others—again reflects capillary 

d viscous effects at low heads. If the latter factors are ignored, and if the 
asic weir equation is written in terms of a discharge coefficient Cg, 


ae: 2 c,V% 3/2 (6) 


: 


/ will be seen that the essential aspects of the Bazin and Rehbock formulas 
ay be reduced for comparison to 


| 2 
j h 

Co =O s pas 
) q = 0-608 + 0.334 (4) (7) 
| 
nd 
h 

C, = 0.605 + 0.08 = (8) 


nich are plotted against h/(h+w) in Fig. 1. The two curves evidently differ 
y little more than 1% for heads as great as twice the height of weir. 
| When Mises,(7) in 1917, determined the coefficient of contraction analy- 
cally for various two-dimensional forms of conduit outlet, he noted that the 
llues for a symmetrical slot in a normal plate at the end of a conduit would— 
Eb combination with Eq. (2)—duplicate the trend of the Rehbock formula very 
sely over a considerable range if the ratio of slot dimension b to conduit 
mension B were used to represent the ratio h/(h+w). His curves for Ce; and 
y are also shown in Fig. 1. To be noted in particular is the limiting the- 
Jetical value 7/(7+2) = 0.611 in comparison with the empirical values 0.608 
id 0.605 chosen by Bazin and Rehbock. One of Mises’ students, Lauck,(8) 
rly a decade later approximated not only the discharge coefficient but also 
Ce nire nappe profile for the weir of infinite relative height, by means of 
Cauchy integral theorem. 


Inspection of the Bazin and the Rehbock formulas will indicate that the lat- 
+ must approach the limit Cg =oo as the relative height of the weir ap- 
aches the limit zero (i.e., the free overfall), a trend that is necessarily 
red by the orifice function of Mises, whereas the Bazin function has the 
lite limit Cg = 0.942. Now it is obvious that discharge over the free over- 
lis finite rather than infinite, since the effective head is the critical depth. 
ver, it was noted by Boss(9) that the critical stage of flow is already es- 

lished upstream from terminal sills well before they attain zero height, 
fich would yield a trend contrary to that of Eq. (6). That is, following the 
cedure of Rouse,(10) the discharge coefficient for such conditions can be 
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determined by introducing into Eq. (6) the critical-depth relationship 


q =-/ gy. = ./g(h+w)? (9) 


Thus, 


(10 


which is plotted at the upper right of Fig. 1. Evidently, the discharge char- 
acteristics of terminal sills follow quite a different functional relationship 
from that of even moderately low weirs, with an intermediate maximum value 
of the discharge coefficient. Now the coefficient of contraction C, can be de- 
termined from simultaneous solution of Eqs. (10), (6), and (2) and plotted as 
at the lower right of Fig. 1. If, after Rouse, this coefficient is assumed to 
approximate the relative thickness of the nappe at its highest point, then this 
plot indicates that the flow profile will also attain a maximum Stage at the 
borderline between weirs and sills. 

In order to investigate these aspects of the problem, the second author 
undertook in 1931-33 in the River Hydraulic Laboratory of the Department of 
Civil and Sanitary Engineering at M.I.T. a systematic series of profile (as 
well as velocity and pressure) measurements on the nappe of a weir that was 
successively reduced in height until the overfall limit was reached. The re- 
sults of the profile tests were reflected in a subsequent book,(11) but they 
were found to cover the transition range inadequately and hence were never ~ 
published in their entirety. Similar measurements over a less extensive } 
range of the parameter h/(h+w) were undertaken later by the Bureau of Re- 
clamation and published(12) for purposes of spillway design. Because it was 
believed that certain characteristics of the weir-sill range were nevertheles: 
of technical interest, generalized plots of the M.I.T. data showing a parallel 
between the discharge- -coefficient and nappe-profile functions were included 
in the Proceedings of the Fourth Hydraulics Conference(13) of the Iowa Insti- 
tute of Hydraulic Research. ; 

Continued efforts have since been made at the Iowa Institute to complete © 
both the analytical and the experimental phases of the weir and sill problem. 
The relaxation method has already been applied to the improvement of 3 
Lauck’s approximation of the nappe profile for the infinitely high weir(14) 4 
at the time of writing the study is being extended to other values of h/(h+w). 
However, Since numerical integration is not as readily applicable to low w 
and sills. as to high, the gaps left in the M.I.T. study were in 1956 filled in 
perimentally by the first author(15) under the direction of the second, to the” 
end of providing a basis for both design and further theoretical analysis. It 
is to the presentation of the composite results that this paper is devoted. 


Experimental Procedure 


The M.I.T. experiments were conducted in a 50-cm glass-walled flume, ’ 
the rate of flow through which was measured with a ventilated sharp-creste 
weir 40 cm high in accordance with Eq. (5). The experimental weir was 
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nitially of identical characteristics, but its effective height was changed by 
he addition of successive false floors on the upstream side. These were 
onstructed of troweled cement mortar over sand, and each consisted of a 6- 
t. level section preceded by a very gradual rise from the original floor of 
he flume. Three weir heights—40, 20, and 10 cm—were produced in this 
vay. The 40-cm weir was then removed; and still another false floor was 
dded at a height of 37.5 cm above the flume floor, with a brass-angle end to 
imulate the free overfall. A supplementary weir crest could be attached to 
e end with wingscrews to produce additional weir heights of 5, 2 1/2, and 
1/4 cm. Full ventilation was provided below the nappe in all instances. 
Profile measurements were made at 1-cm sections upstream and down- 
ream from the crest with point and hook gages reading by vernier to 0.1 
m. For each of the seven weir heights the same nine rates of discharge 
ere investigated; these corresponded to heads on the measuring weir rang- 
g by 2-cm increments from 2 to 18 cm. Supplementary measurements of 
essure distribution along the floor, up the weir plate, and through the nappe, 
S well as checks upon the constancy of the total head at various points, were 
So made, but neither the experimental techniques nor the results are per- 
nent to the present discussion. 
| The tests conducted at the Iowa Institute employed the same type of free- 
werfall structure as the one finally constructed at M.I.T., except that it was 
ilt of brass plate and angles and mounted in a 1-ft. flume. The level por- 
on was 4 ft. in length, its height above the floor of the flume was 1 ft., and 
xe approach curve had the form of a quarter ellipse with an axis ratio “of 
5:1. The adjustable weir plate was machined to a sharp upstream edge, the 
»p surface having a width of about 1/32 in. To minimize capillary distortion 
‘the lower surface of the nappe at small weir heights, the crest surfaces 
re rubbed with paraffine. For the very smallest height (0.005 ft.) the crest 
S remachined to a knife edge. ,The test procedure was essentially the same 
that at M.I.T., although the range of head was limited to 0.1 ft.ch<0.5 ft. 
nd the range of weir height to 0.005 ft<<w<0.05 ft. Discharge measurement 
effected by means of a 90° elbow meter in the 6-in. supply line that had 
: previously calibrated against a 12-in. sharp-crested weir mounted in 
e flume itself, the Rehbock formula being used as reference just as in the 
LT. experiments. 
| 
] 
| 


Discussion of Results 

| 

. As has long been realized, the head on a weir must be measured a con- 
lerable distance upstream from the crest, where the effect of the curvi- 

i ae of flow upon the surface elevation is still negligible. A distance 

ual to three or four times the head is generally considered to suffice. Un- 

unately, with relatively low weirs the effect of channel resistance becomes 

ipreciable, so that the nappe profile is asymptotic to a sloping rather than a 

irizontal line. In other words, the head measurement then is subject to two 

es of error, the one increasing with proximity to the crest and the other 

h distance from it. 

In order to determine that section at which the net error would be a mini- 

m, so that the heads for different weir heights could be evaluated con- 

tently, a composite dimensionless plot was made of all profile measure- 

ents upstream from the crest for the most extreme case—the free overfall— 
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as shown in Fig. 2. Since the computed critical depth was used as the refer- 
ence length, and since the critical section has been shown to be the last at 
which hydrostatic distribution of pressure (i.e., absence of curvilinearity) 
can be assumed to prevail, (10) the abscissa of the desired section should be 
that having an ordinate value of unity, or x/h = 4. Inasmuch as this value for 
the free overfall agrees with that commonly used for weirs of great relative 
height, it was adopted for all subsequent evaluations in the intermediate rangi 
as well. 

From the experimental determinations of discharge and surface profile fo: 
the various heights of weir and sill, the coefficient of discharge was evaluate 
according to Eq. (5) and plotted as in Fig. 3. Herein the Rehbock parameter 
h/w is used as abscissa in the weir range, because of the linearity of the re- 
sulting function. Its reciprocal w/h then logically serves as the abscissa in 
the sill range, the relative scales of the two being so chosen as to yield ap- 
proximate agreement in the intermediate range. All of the Iowa data are 
shown, together with a comparable number from M.I.T. (i.e., all at the highe1 
heads). Inspection of the diagram will indicate that Eq. (8) for weirs applies 
with good approximation over the range 0<h/w<6, and Eq. (9) for sills over 
the range 0<w/h<0.06. Although some discrepancy between the two sets in 
the zone of overlap is apparent, a smooth transition between the two functions 
can readily be drawn. 

In Fig. 4 are superposed a systematic series of nappe profiles extending 
through the weir range to the case of maximum coefficient, h/w = 10. Figure 
5 shows the corresponding family of profiles for sills. Like the coefficient 
plot of Fig. 3, the profiles (except for the lower surface of the weir nappe) 
are seen to be characterized by similar trends, the case h/w = 10 or w/h = 
0.1 evidently representing a maximum for both the discharge coefficient and 
the nappe deflection. 

As had previously been noted by the second author and used as the basis 0 
the generalized profile diagrams already mentioned,(13) plots of z/h as func- 
tions of h/w and w/h with x/h as parameter bear a striking resemblance to © 
the coefficient function shown in Fig. 3. Revised in accordance with the sill 
data now at hand, these diagrams are presented as Figs. 6 and 7, together 
with measured points for two representative values of x/h. The confocal rec 
tilinear portions of the curves below h/w = 5 are identical with those previou 
ly published, but the remainder will be found to differ appreciably. Similarit 
with the coefficient curve is maintained well into the sill range, however, ex- 
cept in the very neighborhood of the overfall limit. 

Although the foregoing experimental data were obtained with heads of veil 
tively small magnitude and hence can be expected to apply only approximate] 
to conditions at larger scale, it is believed that the degree of approximation 
will not be appreciably lower than the experimental precision evident in the - 
figures. The Rehbock equation, for example, indicates that capillary and 
viscous effects change the discharge by about 1% at a head of 0.3 ft. Since — 
this is representative of the heads that were investigated, there is no reaso 
to expect scale effects of a higher order to be present in the results. — k 


CONCLUSIONS 


Through the combination of experimental measurements made at consi 
ably different times and places, information is now at hand for the coeffici 
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of discharge and the form of the nappe profile as the relative height of a 
ventilated, sharp-crested weir at the end of a channel varies over its com- 
plete range from infinity to zero. In the range 0<h/w<6 the coefficient of 
weir discharge is given approximately by the formula 

| 
| h 
| c* 0.61 + 0.08 = 


_ range 0<w/h<0.06, the coefficient of discharge is given approximately 
¥ 


Ca h 


| 
ie 


= 1.06(1 oe 


m the intermediate zone there is a continuous transition between the two 
unctions with a maximum at approximately h/w = 10; the latter might hence 
,¢ regarded as the borderline between weir flow and sill flow. This value of 
/w also marks the condition of maximum nappe deflection, the nappe for both 
veirs and sills gradually rising and thickening as this common value is ap- 
roached. Nappe profiles for design purposes may be determined either by 
terpolation of the profile curves or by direct reading of the generalized 
rofile coordinates presented in the paper. 
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Rehbock 


Fig. 1. Comparison of Discharge Functions for Weirs and Sills. 
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Fig. 3. Experimental Evaluation of Discharge Coefficient for Entire 
Weir-Sill Range. 
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Fig. 5. Variation of Nappe Profile with Relative Height of Sill. 
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Fig. 7. Generalized Coordinates of Upper Nappe Surface for Terminal 


Weirs and Sills. 
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Discussion of 
“FREE OUTLET AND SELF-PRIMING ACTION OF CULVERTS” 


by Wen-Hsiung Li and Calvin C. Patterson 

(Proc. Paper 1009) 
__ WEN-HSIUNG LI, 1aM. ASCE, and CALVIN C. PATTERSON,2 J.M. 
\ASCE—The writers are glad to receive discussion from investigators of 
three hydraulic laboratories. As their comments overlap to a certain extent, 
the authors would like to answer them together. 

The writers use mD for the mean piezometric head at the outlet in Eq. (1) 
and then define an effective mD by Eq. (4). As can be seen from the accom- 
spanying sketch, and from Eq. (3) and (4), 


i) 


: a “3 ~ Vv 
) Effective mD = Actual mD + (A, A.) + (a, a 22g 


ote that the location of the projected hydraulic grade line corresponds to the 
ifective mD, not the actual mD. Thus the values of effective mD presented 
the paper are theoretically the same as those obtained by projecting the 
nydraulic grade line to the outlet section. This can also be seen from 


p 
q. (21) in Mr. French’s discussion by putting ss Gr + Z)dQ = actual mD, 
Lv ; 
| 1 = Qo, and h, + {5-5 2g = Ap. The small differences between A, and A, and 


etween a, and Oo are of no great importance in practice. Since Ay and Oty 


re unknown and will most likely be approximated in practice by Ao and a9 


respectively, the quantity effective mD is even more useful than the actual 


2 

H + sL = effective mD +Ao +o a 
| In presenting the data of effective mD, the writers call the quantity v/V gD 

| Froude Number. Mr. Blaisdell thinks that this quantity is not a Froude 
mber “as ordinarily conceived” relating to free-surface flow. This isa 
atter of definition. It seems that the generally accepted definition of a 
roude number is a dimensionless quantity in the form of v/ VYgD (or its vari- 
ions). It is an indication of the relative magnitudes of the inertia and gravi- 
tional forces(1) in a flow system, e.g. the flow from a closed container 
rough an orifice.(2) Ina particular problem, the Froude number may some- 
mes be interpreted in a particular way, such as the ratio of a characteristic 
locity to the speed of propagation of an infinitesimal surface-wave, or as a 
easure of the ratio of the velocity head to the depth of flow or to the depth of 


| Associate Prof. of Civ. Eng., The Johns Hopkins Univ., Baltimore, Md. 
| Engr., Anderson and Nichols, Consulting Engrs., Baltimore, Md. 
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the conduit. Such interpretation may be helpful in particular cases but, it 
seems to the writers, should not be considered as a definition. The Froude 
number is significant in determining the value of mD because the pressure 
distribution at the outlet section depends on the relative magnitudes of inertia 
and gravitational forces acting on the jet. The writers use v/VgD instead of 
one of its variations to take advantage of the fact that, at the same v/VeD, the 
values of m for unsupported jets are about the same for square and circular 
pipes. 

Considering the nature of the measurements, the differences between 
Mr. Flack’s data and the writers’ data of m for horizontal culverts with no 
outlet-apron are not unexpected. It is believed that the data presented for 
various culvert slopes and outlet conditions are consistent enough for practi- 
cal purposes. 

The writers appreciate the comments on the self-priming action of cul- 
verts. The work of Mr. French in this respect is particularly interesting to 
the writers. They are glad that Mr. French adds Type A to the three types of 
priming mentioned by the writers who, in writing the paper, considered 
Type A as the limiting case of Type 1. It seems that Mr. French’s Type B 
of priming is the transient phenomenon of what the writers call Type 2 prim- 
ing. In the paper, the writers suggest a way of presenting the data on priming 
for practical use. As there is no exact knowledge about turbulent non-uniform 
flow, air entrainment and vortex formation (as of 1957), there is no exact way 
of transferring model data to a prototype. A practicable, though approximate, 
method has to be used. It should be recalled that, in using any empirical in- 
formation, it is a sound engineering practice to allow a margin of safety. The 
data presented in the paper, being experimental results, should be used in the 
same spirit. Friction is a major influence in Type 1 and Type 2 priming. In — 
Type 3 priming where air entrainment is most important and friction only 


1/39 
plays a secondary part, the action depends mainly on D/H. The use of a x 


is not incorrect though not absolutely necessary. 


IO FS ) Ora 


REFERENCES 


1. Engineering Hydraulics, edited by Rouse, Wiley & Sons, N.Y., 1950, p. 140, 


2. Rouse, H. “Fluid Mechanics for Hydraulic Engineers,” McGraw-Hill Co:;m 
N.Y., 1938, p. 19. 


DISCUSSION 1348-5 


Fig. 1. 


ed 


ASCE 1348-7 


Discussion of 
“THE APPLICATION OF SEDIMENT-TRANSPORT MECHANICS 
TO STABLE-CHANNEL DESIGN 


by Emmett M. Laursen 
(Proc. Paper 1034) 


EMMETT M. LAURSEN.!—A discussion such as that of Mr. Shulits can 
only gladden the heart of a writer, for it amplifies and reinforces the original 
contentions set forth. The “formula dilemma,” as Mr. Shulits puts it, is such 
that at the present time the design of stable channels cannot proceed with any 
Substantial degree of confidence on the basis of sediment-transport mechanics 
aione—experience is absolutely essential. Sediment-transport mechanics 
Can serve as a guide and the many formulas can serve as scaling relation- 
ships, however, and in that way supplement experience. When the formulas 
are used to compare different channels the gross quantitative disagreements 
largely disappear and the similarities are emphasized. 

In a personal communication Mr. J. W. Delleur has pointed out that power 
of n in the Brown-Einstein equation is incorrect. In Table I the power should 


be 3 and in Table II the power should be 6/5. 


1. Research Engr., Iowa Inst. of Hydr., Research, Iowa City, Iowa. 
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Discussion of 
“MEASURING STREAMFLOW UNDER ICE CONDITIONS” 


by A. M. Moore 
(Proc. Paper 1162) 


_ CHARLES E. BEHLKE,* A.M. ASCE.—As more and more of the cheaper 
water power sites in the United States are developed, the day approaches 
when American industry must look to the great power potential of northern 
Canada and the Territory of Alaska. In these areas of great water power 
potential the winters may last as long as eight or nine months each year. It 
is in the North that the methods of evaluating ice-affected records of stream 
flow, so well described by Mr. Moore, have greatly magnified significance. 

It is quite common in Arctic and subarctic regions for the winter stream- 
flow to constitute a relatively minor amount of the total annual flow. How- 
ever, Since stream beds may drop quite rapidly, the winter flow can be quite 
important to a hydro development because of the frequent lack of potential 
reservoir volume. It is therefore quite important that streamflow records be 
obtained for the winter months. With a very few exceptions, the methods 
which Mr. Moore has so adequately described of evaluating streamflow from 
Stage records would be quite applicable to these northern areas. 

Discharge measurements by themselves are quite tedious and difficult to 
make in cold weather especially if the thickness of the ice cover on the 
stream is very great. It is quite common for ice thicknesses on northern 
streams to range from Zero to six feet. This presents a considerable prob- 
lem before the current meter can be placed in the stream. The efficiency of 
aman gets quite low as the temperature drops, so a man taking a discharge 
measurement in the very cold weather of the winter may take several days to 
obtain a measurement which would take him only a few hours in the warmer 

. Weather. To add to the difficulty, the daylight hours are much shorter in 
these northern areas than they are further south, so a single discharge mea- 
surement may take two men two or three days to complete if the weather is 
not too cold or otherwise inclement. If the weather becomes colder than 
about minus 40°, very little is lost by simply stopping work and waiting for 
warmer weather. This means that men must wait at a location sometimes for 
long periods of time before they are able to make the discharge measurement. 
It is not uncommon in Alaska and northern Canada for the temperature to drop 

to as low as minus 60°, and colder, which means that during a normal winter 
there is a good deal of time during which the weather is too cold to work out- 
Side efficiently. 

.. _ Since discharge measurements during the winter in northern climates are 

quite expensive, streams may be gaged only a few times during these months. 
The discharge is then known only at the time of the measurement and must be 

“estimated for other times. A record of the stage, even though it is effected 
_by ice cover, anchor ice, or ice jams, would be quite valuable for a deter- 

My mination of the discharge during winter months as has been shown by 


e 
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Mr. Moore’s paper. The obtaining of such a record in a cold climate is, at 
present, practically impossible because of the limitations of the methods 
used for obtaining these records. The standard type of gage well as well as 
the intake pipes leading from the river to the well freeze quite early in the 
fall and remain frozen quite late in the spring and even into the summer. The 
writer has had the experience of being a member of a party which thawed the 
intake pipes to a gage house on the Tanana River, one of Alaska’s major 
rivers, after the first of July. This freezing can affect spring and summer 
records which are very important. In Alaska, since the frost has usually 
penetrated deepest into the ground sometime in May, intake pipes thawed af- 
ter the spring runoff begins may refreeze again very quickly. This points up 
the very serious fact that considerable changes must be made in the matter 
of obtaining stage relationships for these northern streams before good re- 
cords can be obtained. 

The writer would warn the inexperienced engineer against attempting to 
correlate stream discharges with air temperatures at low temperatures. 
Mr. Moore has pointed out that this can be done at moderate temperatures 
below freezing. The writer has seen engineers very carefully adjust their 
estimates of stream discharge with changes in temperature where the tem- 
perature ranges were in the order of from minus 20° to minus 65° with the 
cover on the ground perhaps averaging 3 feet of very light, well insulating 
snow. Since winter stream flow in northern climates is derived entirely 
from groundwater, it is quite improbable that the particular day’s discharge 
from a stream could be affected by the temperature for that day, but it would 
be more a function of a mean temperature extending over a much longer 
period of time as well as the depth and density of snowcover on the basin. 

Some engineers make a practice of installing thermograph thermocouples 
in the gage well instead of in the stream. This can make a considerable dif- 
ference in the temperature which would be recorded for the stream at a par- 
ticular time. The writer has measured differences between gage well and 
stream temperatures of upwards from fifteen degrees F. in the summer. In 
the winter after the gage freezes, the temperature for the stream would 
erroneously be recorded as below freezing by perhaps a considerable amount. 
This is obviously not a true record of stream temperatures, and the use of 
such records should be reviewed critically by the engineer having recourse 
to them. 


STEPONAS KOLUPAILA.*—This paper by Mr. Moore is a valuable contri- 
bution to the complicated problem of determination of streamflow under ice 
conditions. The phenomen of siphon action is of particular interest. The 
writer had occasions to observe similar strange effects in a cross section 
jammed with frazil ice: the water flow, partly in the form of a slurry, was 
concentrated in a few separate channels under extremely different hydraulic © 
conditions. There is no substitute for intelligence and experience in compu-_ 
tations of winter flow, more so than in any other field. 

The writer was very interested in the method of adjustment of daily dis- ~ 
charges applied by the author. Obviously, the author uses the Stout method: 
he corrects the water stages to the “effective heights,” as this is explained, 
e.g., in the “Stream Flow” by Grover and Harrington, p. 285. Another method 
of correction of discharges seems to be of greater flexibility and accuracy: 


* Prof. of Civ. Eng., Univ. of Notre Dame, Notre Dame, Ind. 
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e ratio is determined between the actually measured water discharge Q' 

nd that read on the discharge curve (available for open channel conditions) 

. This ratio k = Q' /Q varies through the winter period due to the changes in 

ickness of ice cover, presence or disappearance of frazil ice, depth of water 

nd, particularly, the smoothing of subsurface roughness by flowing water. 

the obtained ratio is to be interpolated between actual measurements, ad- 

sted at the beginning and the end of ice cover, when direct measurements 

re too dangerous or impossible; meteorological data can be useful in tracing 

€ diagram of k. Daily discharges are to be computed multiplying the value 

ead from discharge curve and the ratio k from diagram. Fig. 1 illustrates 

is method. 

This “discharge correction method” in contrast to the “gage correction 
ethod” was initiated by the writer in 1928 and since that time has been in 

se in many European countries. Several detailed procedures were recom- 
ended for deficient data, especially for previous years without any direct 
inter measurements. Thus, the writer was able to compute the runoff of the 
emunas River in Lithuania for 1812 - 1943, for complete 132 years, one of 

e longest series of existing reliable runoff data. Further developments of 
his method might be mentioned. Mr. A. Ogievskii derived a “hydrometeoro- 
iogical method” for corrections k: air temperature is used; accumulation of 
lhe sub-freezing temperatures aids in construction and detailization of the 
liagram. Recently, Mr. H. A. Klein obtained an approximate empirical formu- 
a for the variation of the k ratio in Elbe River: k = 0.308 + 0.02 T, where T 
s the number of days since ice cover was established. 

This method was presented by the writer at several international confer- 
neces: Leningrad 1928, Tallinn 1928, Seville 1929, Barcelona 1929, 
ielsingfors 1936, Edinbourgh 1936, Lubeck 1938. In this country Mr. N. C. 
qover, chief engineer in the Water Resources Division, U. S. Geological 
jurvey, tried in 1936 to introduce this, as he called it, “Lithuanian method” 
nto practice of his Division. Several valuable reports were submitted on this 
subject, as by M. C. Boyer and B. J. Peterson, G. E. Harbeck, C. C. McDonald 
nd others. Unfortunately, these interesting reports were not published in the 
Vater-Supply Papers. The writer is indebted to Mr. McDonald for a copy of 
is dissertation “Determination of stream flow for periods of ice effect,” 

943, which is of very great value and still retains its importance. 

The writer disagrees with the statement of Mr. Moore that engineering 
iterature contains few references concerning the effect of ice on stream flow 
neasurements. Starting with the pioneer work, done by eminent American 
ydrologists H. K. Barrows, R. E. Horton, and W. G. Hoyt, there are many 
apers in many languages. Some of the most interesting are presented in the 
ubsequent reference list. 

The opinion of the author concerning the accuracy of winter records is 
ompletely true. His suggestion to reexamine all previous winter flow compu- 
ations before 1950 and to improve the annual runoff data is of great signifi- 
ance and should be materialized. The publication of reports on winter flow 
‘ould contribute much to the progress of establishment of the runoff balance. 
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Discussion of 
“TECHNICAL PROBLEMS OF FLOOD INSURANCE” 


by H. A. Foster 
(Proc. Paper 1165) 


G. N. ALEXANDER,* A.M. ASCE.—Mr. Foster’s latest paper concerning 
€ application of statistics to flood estimation entitled “Technical Problems 
n Flood Insurance” is, as he states, based largely on his preceding paper on 
Flood Insurance.”(1) His first article on statistics for engineers, published 

y the ASCE some thirty years ago(2) still remains a classic in this field. 
s the later papers give insufficient technical background to the author’s 
indings concerning flood probabilities it is not possible in this discussion to 
nter into the technical detail that characterised his first article. 

The two basic questions concerning flood probabilities are: 


Question (A) What is the “best estimate” of the probability or return 
period of a flood greater than a specified magnitude? 
To answer this question decisions are required on— 
(1) The type of distribution 
(2) The method of parameter estimation. 


Question (B) What is the “precision” or reliability of this estimate? 


The earlier papers, such as those of Foster(2) and Goodrich(3) dealt only 
rith the first step of Question A, viz. Distribution Type. There is little dis- 
ussion as yet on the second step in engineering journals, but it is treated in 
tatistical text books. 

The importance of Question B was recognised some thirty years ago by 
ne actuary Arne Fisher, in his discussion on Goodrich’s paper\’/) when he 
tated— 


_ “the paper deals only with a subject of graduation, and leaves almost 
‘untouched the far more subtle and difficult aspects of the probable or 
presumptive values of future observations.” 


.S statistics is applied to flood data primarily to enable inferences to be 
lade concerning the magnitude of future floods, this problem is paramount. 
‘his question of precision is closely related to that of sampling errors, which 
s the author states is now being recognised by the engineering hydrologist, 
lthough as yet no particular procedures are generally accepted. 

Given sufficient data and using recognised statistical techniques, Ques- 
ons A & B above can be answered for recurrence intervals not greatly ex- 
eeding the length of record; the difficulties arise when extrapolation is 
ecessary. For an assumed distribution, the precision of an estimate de- 
reases approximately as the logarithm of the recurrence interval and in- 
reases as the square root of the length of record. Consequently extrapola- 
on of short records can give quite misleading results, but with long records 


Engr.-in-Charge of Water Resources Section, State Rivers and Water Sup- 
ply Comm., Melbourne, Australia. 
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some extrapolation is possible without undue loss of accuracy. 

Failure to recognise the magnitude of the sampling errors invariably 
associated with short records was largely responsible for discrediting the 
use of probability methods in relation to floods. The writer believes that 
statistical techniques properly applied must ultimately form the basis of es 
timating flood probabilities which in turn are essential to actuarial estimat 
of premiums for flood insurance. 

We next examine more closely, but still very generally, the means of re. 
ducing the arbitrariness in current procedures in seeking to answer Ques- 
tions A and B. 


Question (A) Methods of Curve-Fitting 

As the author states “the various proposed methods of curve-fitting 
generally produce reasonably good results with a given flood record when 
used to express probabilities or recurrence intervals within the length of 
record.” 


Step (1) Type of Distribution: 

Attempts have been made by some authors to justify certain types of dis 
tribution for flood estimation on “a priori” grounds. Thus, floods are re- 
garded by Gumbel(4) and others as “extreme values” and hence they favour 
the Fisher-Tippett type of distribution. Others contend that the log-normal 
distribution is basically more appropriate, as floods may be regarded as th 
product of a number of independent causes. Arguments of this sort carry 
little weight, as they are seldom backed by study of causes of floods and 
their hydrographs; also the conventional use of the Chi-squared test for 
“goodness of fit” is of doubtful utility. 

Chow has shown(5) that for a specific value of the variance, the asympto 
extreme value distribution and log-normal distribution when standardized a 
practically identical. From unpublished analyses carried out in the writer’ 
office, it appears that the log-normal can also be used in lieu of the finite « 
treme value distribution based on a normal parent by selection of the appre 
priate parent size. This suggests that the log-normal distribution is suffi- 
ciently “flexible” for most cases. The statistical characteristics of the 
log-normal are becoming better known and its conservatism in extrapolatic 
renders it particularly suitable for estimating design floods for spillways. 


Step (2) Estimation of Parameters: 

As long as there are insufficient grounds for differentiating between par 
ticular distributions and methods of parameter estimation, there is some 
justification for the time-saving process of fitting curves “by eye.” Howev 
if extrapolated values are quoted, it is desirable to specify the type of dist 
bution used and the method of fitting, for the different procedures give at 
times very different answers. 

Having selected the distribution, a decision is now required on the most 
appropriate method of estimating parameters. Two methods are commonly 
used by statisticians viz. the “method of moments” and method of “maximvu 
likelihood.” For certain theoretical reasons the maximum likelihood meth 
is preferred by most statisticians, at least for fitting within the range of ot 
servations. However, if curves have to be extrapolated the method of mo- 
ments automatically weights the extreme values, and there may on this ac- 
count be some justification for its use—but this is a contentious technical 
problem outside the scope of this discussion. Even when fitting a log-norn 
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Stribution, estimates from extrapolated curves will differ significantly de- 
nding on whether the fit is to the original observations, or—as appears 
ore desirable—to their logarithms. 


estion (B) Precision of Estimates 

The author refers to the question of the reliability or precision of esti- 
tes towards the end of the paper stating “the question of the reliability of 
robability estimates of future floods is of much greater significance in a 
rogramme of flood insurance than in justifying a flood protection project.” 
y estimates are of much greater value if some indication is given of their 
eliability and this is particularly true of flood estimates. 

This question of the precision of estimates as a function of sample size is 
basic problem of statistical inference. It may be said to date from “Stu- 
ent’s” work on the t-test published nearly 50 years ago. However, there is 
till controversy concerning the more technical problems of statistical in- 
rence such as the relation between the Neyman-Pearson concept of “con- 
idence limits” and Fisher’s “fiducial limits?”(6) In addition to these two 

ethodological developments, “tolerance limits” may also be obtained over 

e€ range of observations. Hydrologists seeking to apply these procedures 
oncerning probability limits naturally have some doubts in trying to decide 
yetween them—apart from any mathematical difficulties that may arise. 
Although the practical differences in the results using alternative methods 
may be small within the range of observations, the validity of inferences over 
he extrapolated estimates is questionable, particularly when the fitted curve 
lepends on the distribution type assumed. For example, when the log-normal 
listribution and the Pearson Type III were both fitted to 50 years of annual 
maxima of monthly discharges of the River Murray at Jingellic, it was found 
hat the log-normal estimate for 100 year return period was very close to the 
ipper 2 1/2% risk curve (95% confidence limits) for the Pearson Type II. 
uch differences, although disturbing to those unfamiliar with the character - 
stics of statistical methods, are believed to be no greater than differences 
.ssociated with other more popular methods of floods estimation for long 
‘ecurrence intervals such as the “maximum possible precipitation” technique. 
The rather misleading title to this technique implies that no particular 
robability is associated with the estimate obtained, nor is its precision in- 
licated. The method does show how, given certain meteorological conditions, 
loods can occur whose recurrence intervals are so great that probability 
nethods which ignore the effect of sampling errors need re-interpretation. 
t is believed that this meteorological method and that of statistical inference 
rom flood data will be reconciled when the probabilities associated with such 
teps as storm maximization and transposition can be ascertained. 


CONCLUSION 


‘It is clear that there is still much to be done both statistically and hydro- 
ogically before the estimation of future floods by statistical inference is 
laced on a firm foundation. With the recent enactment of flood insurance 
sgislation in U.S.A. and the consequent need for actuarial determination of 
remiums, it may be expected that studies of this sort will be accelerated. 
Ir. Foster has made notable contributions in this field, and it is hoped that 
nis paper will stimulate others to publish their findings. 


; 
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STEPONAS KOLUPAILA.*—In his discussion of the determination of floo 
probability, the author states: “After studying and comparing various 
methods of curve-fitting that have been proposed for hydrological studies, tl 
author has not been able to reach any definite conclusions regarding which 
one is the best for use in flood-insurance studies.” As it is known, Mr. 
Foster offered his method in 1924 and prepared his tables, based on two 
characteristic values—coefficient of variation and coefficient of skew. Fos- 
ter’s method was quite popular and widely used here and overseas. After 
some unfortunate failures occurred in extrapolation of probability curve this 
method was almost entirely eliminated. The writer feels that the cause of 
such a failure was the improper application and lack of experience. The co 
efficient of variation can be easily established and checked, although it may 
be based on short period of available data. The coefficient of skew is more 
intricate: the usual period of 20 to 30 years is definitely not sufficient for 
establishment of any reasonable value. Sometimes the addition of just one 
year causes a great change of this coefficient. Therefore it is more sound | 
presume a value by analogy with some more reliable long range data. The 
adjustment suggested by A. Hazen for inadequate length of records is not 
satisfactory. Different rules were offered for improvement of the coefficie! 
of skew. Russian hydrologists made some significant progress: they intro- 
duced different empirical relationships and an additional “guarantee-factor. 
The results obtained with these corrections are more satisfactory. The 
writer has a successful experience with Foster’s method and therefore hope 
that this excellent American method will be restored at home in near future 
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Discussion of 
“FREQUENCY ANALYSIS OF STREAMFLOW DATA” 


| by David K. Todd 
(Proc. Paper 1166) 


CORRECTIONS—Through an oversight, Figs. 10 through 16 were not pub- 
lished with this paper when it was published in February, 1957. These 
figures, therefore, are reproduced herewith. 


M. A. BENSON.!—The author is to be commended for attempting to bring 
before engineers a summary of information on the application of statistical 
analysis to streamflow data. The subject is not entirely new to engineers. 
For a comprehensive discussion of flood-frequency practices and a complete 
bibliography on the subject the writer would like to call attention to the recent 
article of Chow2 on this subject. 

_ Mr. Todd refers to available data in the annual Water-Supply Papers of the 
U. S. Geological Survey. Mention should also be made of the compilations of 
records of surface waters of the United States through September 1950, now 
being published in a series of volumes, each of which covers a major Araane 
age basin of the United States. Volumes already published are in Water- 
Supply Papers 1301, 1311, 1313, 1316, and 1317 respectively. These include 
listings, in readily available form, of annual peak discharges. 

The author states that, “For temporary or small structures where designs 
are based upon flood frequencies about 10 years or less, the basic flood 
method (partial-duration series) should be used;---.” The author repeats 
here the recommendation of the Sub-committee of the A.S.C.E. Joint Division 
Committee on Floods. However, experience has shown that it is preferable in 
most cases to use the annual-flood method, even for small design frequencies. 
The annual-flood method is simpler, the problem of the independence of 
floods is not then involved, and it is always possible to convert from the 
annual-flood to the partial-duration series (see Langbein, reference (22) in 
original article). 

The author says “Most extrapolation by eye of frequencies is arbitrary, 
hence it does not provide a sound basis for estimating frequencies of high 
floods. To overcome this difficulty analytical techniques have been developed 
for fitting a cumulative frequency curve to specified types of frequency dis- 
tributions.” He says further that “With the aid of statistical parameters de- 
termined from the streamflow data and special plotting papers, it is possible 
to estimate say the 10,000-year flood from a record of 20 years duration. 
That such practices are dangerous has been well emphasized by the ASCE 
Subcommittee on Review of Flood Frequency Methods.” 


1. Hydraulic Engr., U. S. G. S., U. S. Dept. of the Interior, Washington 25, 

=D: C. 

2. Ven Te Chow, Hydrologic Studies of Floods in the United States, Symposia 
- Darcy, Dijon 1956-Tome III, Union Geodesique et Geophysique 
ecernationale, p. 134-170. 
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FIG. 12 DAILY, WEEKLY AND ANNUAL-FLOW DURATIO 
CURVES, BIG WALNUT CREEK AT REES, OHI 
(AFTER 10) 
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FIG. 13 FLOW-DURATION CURVES FOR FOUR 
CALIFORNIA STREAMS (AFTER 27) 
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FIG.14 OBSERVED AND FITTED FREQUENCY DIS- 
TRIBUTIONS OF MEAN DAILY FLOWS, 
CUCAMONGA CREEK AT UPLAND, CALIF. 

- (AFTER 27) 


FREQUENCY 


ASCE DISCUSSION 1348-27 


LOG DISCHARGE 


FIG. 15 FREQUENCY DISTRIBUTION OF MEAN 
DAILY FLOWS, FISH CREEK NEAR 
DUARTE, CALIFORNIA. (AFTER 27) 
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FIG. 16 FREQUENCY DISTRIBUTION OF MEAN 
DAILY FLOWS, SALMON RIVER AT 
SOMESBAR, CALIFORNIA. (AFTER 27) 
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The writer questions the implied superiority of fitted mathematical dis- 
tributions over graphical techniques. It has not yet been possible to demon- 
strate a priori that flood peaks should conform to any mathematical or sta- 
tistical distribution. Hence, it is not only dangerous but meaningless to 
estimate a 10,000-year flood from a 20-year record by means of a fitted 
curve. Pertinent statements on the general subject of graphical versus 
mathematically-fitted curves have been made by Ezekiel:3 


“When there is some logical basis for the selection of a particular 
equation, the equation and the corresponding curve may provide a definite 
logical measurement of the nature of the relationship. When no such logi- 
cal basis can be developed, a curve fitted by a definite equation yields on- 
ly an empirical statement of the relationship and may fail to show the true 
relation. In such cases a curve fitted freehand by graphic methods, and 
conforming to logical limitations on its shape, may be even more valuable 
as a description of the facts of the relationship than a definite equation and 
corresponding curve selected empirically. 

“In any event, estimates of the probable value of the dependent variable 
cannot be made with any degree of accuracy for values of the independent 
variable beyond the limits of the cases observed; and can be made most 
accurately only within the range where a considerable number of observa- 
tions is available. It may be possible to extrapolate the curve if its equa- 
tion is based on a logical analysis of the relation as well as on the cases 
observed; but in that case the logical analysis, and not the statistical ex- 
amination, must bear the responsibility for the validity of the procedure.” 


See also Thomas,4 who says, “Indeed insofar as extrapolation made be- 
yond the range of the recorded floods is concerned, the most elaborate analy- 
tical procedure of curve-fitting gives results that are no more reliable than 
these obtained by a simple estimation by eye of the flood frequency curve on 
any kind of probability paper.” 

The author refers briefly to the method of regional flood-frequency analy- 
Sis which can be employed to estimate flood frequencies in ungaged areas. 
His figure 6 illustrates such methods, developed for the Atlanta, Ga. area, but 
he omits mention of the many statewide regional-frequency studies in which 
such methods have been described and developed for a large part of the 


- country. 


The author discusses the fact that in some regions floods may be caused 
by either rain or snow melt and, because “two different populations are inter- 
mingled----- it follows, therefore, that not one but two straight lines of differ- 
ing slope should be used to fit such flood frequency curves.” Because the 
populations are in truth intermingled, the use of two straight lines does not 
follow as a solution. There is in general no particular flood magnitude below 
which all floods are caused by one factor and above which they are caused by 
another. In most parts of the country floods have multiple causes which are 
not distinct but combine or merge one into the other. The situation has been 
well analyzed by Hazen,° who many years ago analyzed the effect on 


3. Mordecai Ezekiel, Methods of Correlation Analysis, 2nd Edition, 1941, 


John Wiley & Sons, N. Y., p. 127. 
4. Harold A. Thomas, Jr., Frequency of Minor Floods, Journal, Boston 


Society of Civil Engineers, October 1948, p. 428. 


5. Allen Hazen, Flood Flows, John Wiley & Sons, N. Y., 1930, p. 102. 
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frequency curves of the multiple causes of floods, but said, “It is true of 
course that these causes grade into each other and that no lines of demarca- 
tion can be drawn.” a 
The discussions of the frequency analysis of low flows and of allflows 
could be expanded considerably, i a survey of these subjects is to be at all ‘ 
complete. 7 
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Discussion of 
“BUTTERFLY VALVE FLOW CHARACTERISTICS” 


M. B. McPherson, H. S. Strausser and J. C. Williams, Jr. 
(Proc. Paper 1167) 


TURGUT SARPKAYA, | J.M. ASCE.—Valuable information adding to the 
gradualiy increasing knowledge of the performance of butterfly valves is con- 
tributed by the authors. The paper constitutes a well-arranged review of 
previous studies in addition to the newly presented basic analytical relation- 
ships, their interpretations and sensible applications to actual flow condi- 
tions. However, the authors have not dealt with the following features of the 
problem: (a) Variation of contraction coefficient at each end of the valve 
plate with the angle of valve opening and with the angle of complete closure; 


_(b) Stability of flow for various valve openings; and (c) An explanation for the 


large scatter of experimental data for small values of o& (@<20°) or for near- 
full-open positions of valve plate. 

For some time, this writer has been interested in the application of mathe- 
matical methods to practical problems of fluid flow, and has contended that 


'many existing semi-empirical relationships would find wider application if 


all parts of the relationships were, wherever possible, freed from empirical 
values and experimental coefficients. Having once obtained the pure mathe- 
matical results, engineers can then compare them with the experimental re- 


sults. Since the writer feels that the results obtained by the authors are sig- 


hificant, it seems worthwhile to apply the methods of hydrodynamics to 


butterfly valve flow in order to have a better understanding of the variations 


of various parameters involved in the phenomenon. From this viewpoint, the 
omission of any discussion of the rather surprising analogy which permits the 
use of hydrodynamic theory for certain problems of viscous motion is deemed 
to be regrettable. Having this object in mind the theoretical solution of the 
present problem is given in this discussion. Then the experimental results 
are compared with the theoretical ones and answers provided for the points 
Stated above. 

It is well known that many plane irrotational flow patterns which include 
the formation or deflection of free jets can be analyzed completely by the 
method known as the Helmholtz(1)-Kirchhoff(2) theory of free streamlines. 
The results obtained may be significant not only for the corresponding two- 
dimensional flows of real incompressible fluids but, in some instances, for 
their three-dimensional counterparts as well. The analyses consists of the 
definition of successive conformal transformations involving a hodograph, or 
velocity plane, and the application of Schwartz(3)-Christoffel(4) transforma- 
tion. The writer would like to emphasize the fact that from the mathematical 
viewpoint, the direct calculation of free-streamline flows is restricted to two- 
dimensional irrotational flows of incompressible fluids which are free from 
gravitational effects and for which the solid boundaries are straight. As flows 
of real fluids seldom satisfy fully these stringent restrictions, the practical 
value of the results is open to question. The significance of the various 


i. Nazilli, Turkey. 
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restrictions can be shown only by comparing the calculated results with those 
obtained from observations of real flows. For the present problem since the 
shear is not a primary factor, (the effect of shear will be discussed later) 
for the flow accelerates rapidly, the differences between comparable ideal 
and real flows are small. The importance of gravitational effects for the 
coefficient of contraction is already known to be definitely secondary. Being 
fully aware of the restrictions and limitations of the applicability of this 
method, the solution of butterfly valve problem proceeds as follows: If the 
coefficients of contraction at points B and H are Cc 1 and Cc2 respectively as 
is shown in the physical plane (Fig. 1), then the main objective of the present 
investigation becomes the determination of the relationships: 


Saf, ° 


The existing relationships are: 


n=, 4 (13g ) 

m= Git (I-sne 7 

fast, =a (4a), (4b) 
one 


(5) 


a 
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in which n-and m-represent the thickness of the jets at B-and H-respective- 

ly; 7 is the width of channel; B the angle of closure; o the variable angle; and 

L the length of valve plate. The velocities of both jets are equal to Vj and the 

velocity far upstream from the valve plate is Vo: 

If z-and w-are used for the complex variable and the complex potential, 
then the variable ¢ defined as (dw/dz) is given by 


= (6) 
$=-u+iva—qe 


in which u-and v-are the velocity components in the x-and y-directions, q is 
the magnitude of the velocity vector, and @ is the angle of its inclination. 
If the variable Q, 


V; 
Q=In ae =-lnj rie (7) 


is introduced, the flow region in the physical plane (Fig. 1) can be mapped on- 
to a corresponding region in the 2 -plane, as indicated in Fig. 2. Either ¢ -or 


bette 


z HHH oe me Li 


Q-must be related to w by means of known methods through still another com- 
plex variable designated as t. This relationship completes the rather com- 
plicated functional analysis because, as has already been seen, either ¢-or 
Q-can be expressed in terms of (dw/dz). P 
Any polygonal boundary in one complex plane can be transformed into the 
entire real axis of another through application of the Schwartz -Christoffel 
theorem, the interior of the polygon being transformed into either the upper 
or the lower half of the second plane. Thus, the outline in the -plane is 
directly transformable into the real axis of the t-plane. Therefore, ut craig 
out the 2-plane and lay it along the real axis of the t-plane as shown in Fig. 3, 
the interior of the Q-plane transforms into the lower half of the t-plane. The 


transforming function can be written as follows: 
‘ 


q 
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o= "Trt Ga 
-oo -f -l -a +I d +> 
t. Plane 


FIG. 3 


As many factors (t - ty) are introduced into the transformation as there are 
vertices of the polygonal boundary in the 2-plane. The value tp is the locus 
of the transformed point on the real axis in the t-plane, and 7, is the exterior 
angle of the polygon at that point in the t-plane. Three of the tn values can be 
assigned arbitrarily but with due regard to the geometrical characteristics of 
the Q2-plane, and the remainder must be given parametric values to be evalu- 
ated in terms of given quantities, along with M-and N, from the integrated 
function. 

The desired mapping function between 9-and t-becomes: 


uf (t+a) dt +N (9) 


M 
V+) (t41) (4-1) (td) 


The transformation equation between the w-plane and the t-plane is ob- 
tained directly from the complex potentials for sources and sinks, point A be- 
ing considered as a source of strength (27)V)) and points B and H as sinks of 
strengths (2nVj) and (2mVj) respectively: 


= ¥ {nin (t-1) +m ln (t+!) — los y ln(t+a)} (10) 


or 


wali fs fe mp ssh -—(Ga G,)ln(t+a)} (10a) 


sas} 


and w-is connected to z by the relation, 


we olde hee (11). 
z= ‘fe dw E 


The evaluation of this and the other foregoing integrals are shown in the 
appendix of this discussion paper in order not to lose sight of the practical 
nature of the problem under consideration. 

The results or the functional relationships obtained from these equations 
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aa 
a8 


C rare & f Tet 
are shown in Fig. 4. The significance of this figure is obvious. The contrac- 
tion coefficient C,j varies greatly witha, but very little with the angle of 
closure 8, whereas the contraction coefficient C.2g varies appreciably with 6, 
put very little witha, for the values of a greater than approximately 159°. 

But for small values of a, Cc2 varies extremely fast and approaches unity. 
The values of the corresponding angles (corresponding to@<10°) used in 
valuating the elliptic functions (see appendix) approach their limit so closely 
hat difficulties are encountered because of limitations in the tables available 
0 the writer.(5,6,7,8,9,10) For those portions of the curves which could not 


7 


ese 
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be computed readily the trend is indicated by broken lines. 
In accordance with the authors’ equations (1) and (7a), the discharge coef- 
ficient CQ for free discharge can be written as follows: 


(12) 
C= FE (Git Con) Se 


C@ values computed from Eq. (12) for 8 = 75° using the contraction coeffi- 
cients obtained in Fig. 4, are plotted in Fig. 5, and the solid-line is obtained. 
In the same figure the experimental data obtained by the authors for 4"'CDC 
valve and 6"CDC valve are also plotted. Though, the complete closure angle 
of the latter valve differs 2.59 from the former, another theoretical curve has 
not been drawn for 8 = 77.5° in order not to make the figure illegible. The 
agreement between the theoretical and experimental results is fairly good, 
the theoretical discharge coefficient being somewhat smaller. 

A similar computation has been made for 8 = 909, and the resultant the- 
oretical curve with the experimental points (Genissiat and Dow data, as re- 
ported by the authors) are shown in Fig. 6. The agreement between the the- 
oretical and experimental results is excellent. 

So far we have discussed the effects of geometric variables on each coef- 
ficient of contraction and on the discharge coefficient. It is seen that the 
coefficients of contraction are functions of @-and B-only. Under the heading, 
“General Analysis for Enclosed Flow,” authors state that: “In these tests, 
the magnitude of y (see authors’ Fig. 1) was not noticeably changed by varying 
the discharge.” Authors are perfectly correct in their observations, because, 
theoretically the value of y (n-and m-in Fig. 1 of this discussion paper) does 
not depend on the magnitude of discharge. 

The second point we have proposed to discuss was the stability of flow for 
various valve openings and in connection with this the explanation for large 
scatter of experimental data for very small values of w-or for near-full-open 
positions of valve plate. The stability of flow depends on the stability of the 
stagnation point on the valve plate. As @-becomes smaller, the stagnation 
point or rather specifically speaking the stagnation line approaches the lower 
edge of the valve plate and theoretically it reaches the edge when a-is zero. 
But, practically the stagnation point or line reaches the sharp end of the plate 
when a@- becomes about 15°. To illustrate this, a curve has been drawn com- 
puting the location of the stagnation line for B = 90° and for various values of 
a, as shown in Fig. 7. The dimensionless parameter A -becomes practically 
zero for a = 20°. 

In the case of axially symmetrical flow, the locus of stagnation points is 
curvilinear. Therefore, the stagnation points in the vicinity of valve axis are 
very close to the edge of valve plate, whereas the stagnation points in the 
vicinity of the diameter perpendicular to the valve axis are relatively far 
from the edge. The writer believes that in the case of axially symmetrical 
flow for a critical value of a -the stagnation points relatively near the plate 
edge jump abruptly to the edge of the plate changing the discharge coefficient. 
In the two-dimensional case the entire stagnation line jumps abruptly to the — 
edge of the plate for small values of a. For a critical value of a, however, — 
the stagnation line is unstable and jump to the edge and then back alternately 
resulting in unstable flow. 

It would be very interesting to compare the discharge and contraction 
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coefficients obtained from two separate tests as follows: (a) First setting 
‘& = 250 and decreasing it a few degrees each time, (b) First setting a = 5° 
dd increasing it a few degrees each time. Each test might well result ina 

eparate curve in the range of critical values of @. The results obtained 

rom such an experiment should prove useful in obtaining a better understand- 
ing of the stability of flow. 
| It seems to be evident from the foregoing explanation that the scatter of the 
experimental data for small values of a-is largely due to the rapid and abrupt 
ariations in the location of stagnation line in that region. Moreover, the 
eectsy of valve plate also has an important effect on the variations of afore- 
mentioned coefficients. Theoretically, the plate is a thicknessless or one 
jimensional surface, but practically it has a thickness. Therefore, as noted 
oy the authors under the heading: “Comparison of Actual and Analytical Flow 
~oefficients,” “For a blade with a relatively wide tip, the point of separation 
vill be at one or the other edge, depending upon the relative blade angle. The 
sudden, slight shift of Co trends for some valves is due to this transfer of the 
separation point.” 

The writer would like finally to comment on the discharge and on the con- 
raction coefficients for submerged flows. It is evident from the experimental 
ta given by the authors that the discharge coefficients for submerged flows 
ire larger than the ones for free discharge as amply explained by the authors. 

ia have observed that as stated under the heading: “Free discharge and 
ubmerged flow characteristics,” “The free-discharge and submerged char- 
— appear to be related, but are distinctly different at near-open 
Ositions. In submerged flow, the fluid fills or nearly fills the void between 
e two jets up to the blade, resulting in an expansion of the given fluid within 
e€ same fluid. In free-discharge of a liquid, the void is normally filled with 
lir. For this reason, it would be logical to assume that the flow characteris- 
ics might not be identical even though the installations are quite similar.” 
The authors, therefore, assume that the action of turbulent shear in 


y 


“2 Tins 
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causing the submerged jets to expand has some influence upon the pressure 
field that controls the rate of efflux. This writer was inclined to accept this 
explanation if it were not due to another contradictory observation and state- 
ment made by Mr. Curtis(11) and his associates on the “Coefficient of con- 
traction for a submerged jet.” Mr. Curtis states that “The minimum dimen- 
sion of the submerged-jet profile is somewhat larger than that of free-jet, 
the coefficient of contraction being correspondingly higher—i.e., Ce = (0.792) 
= 0.63.” At this point one would think that since the discharge coefficient de- 
pends upon the sum of contraction coefficients, the discharge coefficient also 
would be larger substantiating the experimental data given by Mr. McPherso1 
and his associates for submerged discharge case. But Mr. Curtis continues 
to say that “Since the discharge coefficients of the two jets have already been 
seen to be essentially the same, (experimentally), one is forced to conclude 
that the rate of flow through the orifice is essentially independent of the tan- 
gential forces exerted upon the jet by the fluid into which it emerges.” This 
writer is forced to conclude that the departures of the characteristics of sub- 
merged flow case from free-discharge case requires further study and a bet: 
ter explanation. Unfortunately, the authors’ explanations can not be advocate 
as a panacea. 
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APPENDIX 


The integral given by equation (9) can be decomposed into two component 
| integrals: 


t dt a at 
pence. ere 
| V(t7=1)(t-a)(t42) V(t7=1)(t-a)(t42) 
| Let us consider the transformation of Legendre:(12) 
| 2 
’ a sin“p + b, (14) 


2 
cosin’g + 4, 


Ag, by, Co and dg can be so chosen that the above integrals simplify and re- 
duce to known integrals. It should be noted that these integrals, depending on 
the values of t, require separate treatment for each interval under con- 
sideration. 

I—In the intervals DE, and FEg, i.e., fort >d andt < -f, one has 


| 25 
| pits (Sat etna psa OC 182) (15) 


(ae£)sin®g,- (14f) 


4/ (t-4) (142) 


| $,= sin » o€p.£¥2, (16) 
(t-1)(d4+f) 
15 


2(d-1)(d+f)(14f) sin $4-cos #4 
at = $$$ $$$ 7 (17) 
[(ast)sin?¢,- (1+2)] 


If Eqs. (15) and (17) are replaced in Eq. (13), one obtains 
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2. 2M (def )sin?g,- d(1+2) ag, z 
Hist) (19a) (defysin?g,- (142) /1-k2sin°g, 


au d 
teas ee” “ 
(1¢f) (144) 1-ky sin #4 


Integral (18) can be separated into a first and third kind elliptic integral: 


2M 


Qu aon | (448) F(D,, 2/5 ester } (19) 
msl sa) (Pp ep CATS 5 ee 


7 9(9,n,k), is the third kind elliptic integral of Legendre, and may be written 
as 


Tit, 92,%) = I ere aS 


(20) 
(14nsin?9)//1-k*sin¢ 


The parameter n-is taken as real, and may take any value, positive or nega- 
tive. For the intervals under consideration, 


a+f 
py ER —orXny 1 
Legendre’s integral is connected with Jacobi’s by the relation, 


INP See oe eee ee ee 1 CREA SS , (21; 
(#4 sMky) = uy on(e; si, ) an(Ep Ag) eA saiean ae 


in which uy is written for F(D1,,). sn, cn, dn, are elliptic functions and K’ 
is the complete elliptic integral of the first kind for the complemantary modu 
lus. €y is given by the relation 


d+f 2 2 ° 1 
= —— =- ki sn“(€,4+ 1K’, ki) = - 
1 1+f 1 1 1 are (é, “a ay 
and 
2 14+f 2 d-1 2 ad-1 
a+t ioe te Stee aed 


In order to make the notations simpler, sn(€ 1,k;), cn(€4,kj), and dn(e 1,1) 
will be denoted by sn €j, cn€1, and dn €j, respectively. 
The solution of Jacobi’s integral is given by 
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. U, CnNé, dn - 
sné, 2 H(€,+u) 


in which H(€4 ¢ uy) is Jacobi’s Theta-function,(13) and 
2 2(f+a) 
(142) (tea) 
ence Legendre’s integral becomes: 
sné 1 H(E,-u, ) 
Th 6,204 9k,) = - ——— fy 206 94) + —In — 1+ (23) 
. eng, -dné, 2 H(éE, +u,) 


"rom Eqs. (19) and (23), we have 


| 
{ g, = 
| ae {(rearay=ta-n) - [», 2(€,)+ Lin pod ase “Ihe 


| 


hich is the solution sought for the intervals under consideration. The con- 
ants M and N are determined from the comparison of 2-and t-planes for the 
rresponding values and were found to be one and ia resnectively. 

| Il—For the interval HB, i.e., Itl<1, 

tis taken as 


2 rs 
) £@ sa Sal ied Be (25) 


2 sin°g,- (14f) 


d the foregoing manipulations repeated, one obtains 


2M 


V(14a)(142) 


which 


{(a-£78(0%) + (£-1) Th(go9- ’ «,) 26) 


Ne =- ee which varies between zero and axe 
2 1+f 2 


22 
iting Uy = F(g,,k,) and ny = - k,sn (G 9°Ko) =- 
has 


Lt 


0(E,-u,) 
Z(En sk pe 27 
{ (Eq9KQ)+ asta) 2) 
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and therefore, 


ee es ee ee =| 2(é,)+-tn 
es. ( ares (a ke ats ae be (Q)+> 2 “eral 


in which 


fm tnint (142) (14t) 
2 


mn 2(tef) 


ko = 2(f+d) and Ks, = kj = d=1 f-1 
1 Vi (442)(144) (a41) (£41) 
Tse £-1 f-1 
eee PCE 2 Vesa ’ 2 Vet 


6 (E tu) also is Jacobi’s Theta-function and connected to H(e tu) by the relatic 


@(z+ik’) = iBH(z) 
where 


ImB = EK. _ wiz 
4k 2K 


K and K’ are complete elliptic integrals of the first kind connected by the fol 
lowing relationships: 


K = sn7'(1,k) Pa) sn7'(1,k’) = en7*(1, {1-k*) 


I—For the interval BD, i.e., 1<t<d, 
If t is taken as 


2 
(d-1)sin"g, + (1+) (2 


(d-1)sin"g, - (144) 


one obtains 


S2= 


(a-1)uz +2][(¢.,- iris) a 
ese! srcdits E's} AN . 


Repeating the intermediate manipulations, one has 
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: at Picks uae (31) 
ls, (€5)4+ Lin | +N 


Qa atone" yU+ 2 
cnéz.dn&, 6(€,4u) 


v(142)(1+4) 
in which 
? id 
ee me MSRP EST) Cy? ek eke ee 


3 (ast)(t41) 2 3 (142) (148) 


| a ee d-1 2 
| Pomme s. OO ((S5iits iin = 9 75 €ti5 C0 


| $, = sin~? (@)-1) 
| 3 V (a-1)(t41) 


IV—For the interval FH, i.e., -1 >t 
‘tis taken as 


WV 


=f 


2 
d(f-1)sin - f(d+i 
zr guictlh )s By (d+1) 


(32) 
nag 
| (f£-1)sin“g, + (d+1) 
t obtains 
Wiser) (144) : 


« nace » ngs-k,on°(4€,,k,) = kgsn*( €4k¢ )/on”(€, 5X4) 


n,=k7tan“@ » sind = sn(é,,k4) 
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AY x’ @(u,-1€, ) 
Tle u,+ se eee EOC Soe (34) 
dn ( E,2k4) O(u,+1€, ) 


From Eqs. (33) and (34), one has 


(€,k4).cn(E, sk; ) 
Qe tt —ftantng - (abS) eee ee 


(14a)(14f2) an( €, 9k‘) 
ies (35 


Q(u,-LE 
[iu,2(2e,) +i tip ee) 
) 


O(u,+ie, 
in which 
ae Pe 2(d4f) 
u, = F(g,,k,) kf a ko = 
Nhe Sees ties Paleo! (142) (144) 
2 _ (d-1)(f-1) : +1 
ki = sn(&,,k4,) = /— 
4 (as1)(t41) 44 i , 
, a-1 F ft 
en(€,,k,) =/—— dn( é,,k’) =/—— 
494 fad ? ‘49 4) an 
and 


4 


Eo { sinc tan”! 4 | 


Intermediate steps in the solutions of the above integrals have been omitt 
for the sake of brevity. However, the reader will find an excellent account o 
the theory of Jacobian elliptic functions in reference (14), and in the other 
references stated previously. 

In order to obtain z from Eq. (11), dw should be written as, 


Ae sina,, © ¢ Corte 
dw = opi (1- Ol oe Set Olan Oe (3 
af SPI Fe eh pereg er 2 eet , 


The rest of the integration could only be performed by numerical methods. | 
The best and well-known method for this purpose is King’s(15) A.G.M. scale 
method. 

If, for a fixed value of 8, numerical values are assigned to a,f, andd 
satisfying the boundary conditions in the transformed planes, corresponding 
simultaneous values may be obtained for Ccj and Cc2 anda@. It should be 


a 


= 


u 
H 
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noted that this process was a tedious trial-and-error solution. For this rea- 
son, the plot in Fig. (4) was prepared to make possible the intermediate de- 
termination of an approximate result through the practical range of values. 
For large values ofa, say greater than 15°, and for a given B, Cce2 values 
are very nearly equal. Thus only three values of C¢9 for a given B were 
computed and a curve drawn through them. Similarly for a fixed @ and for a 
wide range of values of 8, Ccj values are very nearly equal. Therefore only 
four values of C, 1 were sufficient to draw the appropriate curve. 

The writer believes that a more accurate numerical computation could be 
carried out. However, this has been precluded by the limited time available 
for the discussion. Besides, the above results are deemed to be sufficiently 
accurate for practical purposes. 


G. DUGAN J OHNSON, 4 M. ASCE.—The authors do not concern themselves 
with the characteristics of, or problems associated with, butterfly valves for 
“wide open or fully closed” service in penstocks and at hydraulic turbine in- 
leis, but limit their discussion “to the flow of water through control-type 
valves for near-closed to near-open positions.” They say that “for most ap- 
plications flow characteristics are a function of geometry only” and attempt 
to demonstrate “the effects of geometric variables and present a basis for 
the analytical determination of flow coefficients” from which they indicate a 
method for predicting the inception of cavitation under regulating conditions. 

The authors state that “flow coefficients for the near-open positions are 
‘influenced by the relative blade thickness and the blade shape,” “viscous ef- 
fects are noted witha, = 0° to 20°” and “attainment of a high degree of test 
accuracy within this range is difficult because of the relatively large dis- 
charge rates and small head losses involved.” They also note that discrepan- 
cies in various published data for the near-closed positions “can be explained 
by the..... difference in closure angle.” 

For the intermediate valve disc positions, they developed a two- 


| dimensional plastic flow section (8 wide by 1" thick) to observe and mea- 
| Sure the shapes of the two separate jets issuing from the orifices formed be- 


tween the side walls and the centerline profiles of the discs modeled. They 
reasoned that “since the jet area is a function of the variable (y), measure- 


| ments of the jet thickness in a two-dimensional flow model, representing a 


narrow section of the center part of the blade, should provide a basis to un- 
derstand and predict prototype performance.” The writer agrees that such 

Studies should be of qualitative assistance in understanding the characteris- 
tics of flow patterns at the centerline, but cannot accept the quantitative re- 
Sults as being even approximately applicable to prototype performance. 

It is a well-known fact that the discharge and contraction coefficients of an 
orifice of any shape are critically affected by the edge condition, which must 
be accurately reproduced in the model in order to predict prototype per- 
formance. The authors recognized this when they added the dotted line to the 
right hand orifice in Figure 1, but they minimized it. unduly when they as- 
sumed that centerline flow patterns were even remotely similar to flow lines 
around any circular disc nearer its trunnions. 

_ Consequently, the writer feels that any attempt to predict prototype dis- 
charge performance for valves in pipelines on the basis of such two- 


dimensional model tests would only be adequate where extremely rough 


ii. Chf. Hydr. Engr., S. Morgan Smith Co., York, Pa. 
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approximations were acceptable, and that prediction of cavitation inception 
PRP at: which depends upon the maximum (rather than the average) velocity 
Squared’... 2% on the basis of such tests would be meaningless. 


HENRY VOLTMANN.!—The authors are to be commended for the very in- 
teresting, complete and valuable paper they have presented. Definite knowl- 
edge of butterfly valve performance becomes of vastly increasing importance 
in the hydraulic field and while many reports have been available of past 
tests, the data, as referred to by the authors, have not always been specific 
as to the test conditions. 

It may be interesting to add to the data presented the following information 
concerning a 6" butterfly valve. 

Independent water flow tests were run on this butterfly valve in a “continu- 
ous pipe” installation by Mr. J. J. Colleville under the direction of Mr. M. B. 
McPherson at Lehigh University in 1953. The data from these tests agree 
very favorably with the data in Figure 7 and Table IV of the authors’ paper 
(nomenclature is identical): 


_ hades sos, Wie yy Bangeot Np No. of Test Runs 
10° 1.48 3.3 - 5.9 x 10° 5 

20° 1.00 2.5 - 5.9 x 10° 6 

30° 0.613 2.6 - 5.6 x 10° 6 

40° 0.38 2.5- 5.6 x 10° 6 

450 0.31 ay tee ey ae cao 5 

50° 0.23 0.5 - 4.0 x 10° 11 

60° 0.13 1.3 - 2.7 x 109 7 

70° 0.040 0.5 - 1.0 x 10° 5 

80° (Closure) --- A 


The valve was lens-shaped with a t ,,3,/D of 0.188. There was no indica- 
tion of any variation in the flow coefficient due to viscous effects. The C 
values obtained for individual test runs varied a maximum of 5% from the 
average values given in the above table, except at 700, where only small rates 
of flow were attainable and the errors were consequently higher. 


The following piezometer station—distances were employed in the tests 
(four taps/station): 


Distance, in pipe diameters, from valve body 
Upstream 2D, 4D, 10D, 20D and 30D 


Downstream 1D, 2D, 3D, 4D, 5D, 7D, 9D, 11D, 13D and 15D. 


1. W. S. Rockwell Co., Fairfield, Conn. 
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A straight piezometric head line was obtained for all tests with the up- 
stream stations, to and including the 2D station. The piezometric head line 
was Straight in the downstream pipe, starting at the 5D station for blade 
angles of 10° to 40°, and starting at the 4D station for 459 to 70°. The valve 
loss, AH, was obtained by graphically extending the upstream and down- 
stream piezometric head lines to the valve and measuring the difference. As 
with the authors’ test data, the downstream pipe was filled with water in all 
test runs. 
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RONALD E. NECE.! J.M. ASCE.—The authors have presented the results 
of an interesting and comprehensive series of tests performed upon butterfly 
valves, and the information presented adds to the knowledge of this type of 
valve. The authors have also performed a service in presenting the results 
| of some other investigators in a systematic form. Other data were undoubted- 
| ly considered, but not presented, in this paper. Ee seace ty the results of 
Weisbach, originally given in 1845 and tabulated in (D- 1),2 were the earliest 
given for butterfly valves. 

Of the many types of valves used to regulate flows for different applica- 
tions in continuous pipelines, the butterfly seems to lend itself most readily 
| to a theoretical analysis which will predict, with reasonable accuracy, the 
non-recoverable head loss through the valve. This discussion shall be con- 
| fined to this one area of the investigations presented by the authors. 

The head loss caused by the presence of a constriction in a continuous pipe 
| which causes a contraction of the flow, followed in turn by the rapid expansion 
| of the stream to the full conduit area downstream from the valve, may be 

| analyzed by the combination of energy and momentum considerations. 

This head loss may be expressed, using the authors’ notation, 


Ahe (ve-vs)” (1) 
29g 


| which is the familiar expression for the Borda-Carnot loss, where only the 
| flow expansion downstream from the constriction is considered as being a 


source of loss. 
This expression may be rewritten as 


-(A__ |) _v_ (2) 
an ae I 


Asst. Prof. of Hydraulics, Dept. of Civ. and San. Eng., Massachusetts Inst. 


‘of Technology, Cambridge, Mass. 
2. Numerals in parentheses refer to corresponding items in the list of 


References in this discussion. 


ra 
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where: C, = coefficient of contraction 
a = transverse projection of the restricted area of the valve. 
A = cross-sectional area of the valve body downstream from the 
constriction, in many cases identical with the area of the pipe. 
V = average flow velocity in the pipe. 


Equation (2) may be written in the common form 


ra 
Ane Ki 55 (3) 


where Ky = dimensionless head loss coefficient. For discharges generally 
associated with butterfly valves the Reynolds numbers are sufficiently high 
so that the flow geometry is essentially constant with varying discharges for 
a given valve configuration, with the result that Kj is also a constant. Any 
theoretical expression for Ky, thus depends upon the value of Ce associated 
with the area ratio, a/A. 

The following simple “model,” shown in Fig. 1, may be considered as one 
method of use in predicting head loss characteristics. The valve is con- 
sidered to be in a two-dimensional duct of height B; the disk, which is pivoted 
at the duct center line, is of zero thickness, has the same width as the duct, 
and affects full closure of the valve at an angle @ of 90°. The flow past the 
disk is then divided into two jets of widths Cc, b and Cegab. 

For various angles o@ and slot height ratios 2b/B von Mises has calculated 
the values of Cc for irrotational flow. These are listed in (D-2). From these 
calculated results the average value of C,, for the valve geometry may be ob- 
tained by interpolation for the assumed condition of equal pressures, and 
velocities, at the vena contracta of each jet. The distribution of the flow past 
the two ends of the disk is also determined as a part of this calculation. Thus, 
for the “model,” Ky, may be expressed as a function of o@ alone. 

Limitations to the validity of the “model” are immediately obvious. The 
disks of actual valves possess a finite thickness of both blade and hub; thus, 
when the valve is fully open, there is still a reduction in flow area through the 
valve, with the consequent head loss. Actual disks are usually tapered toward 
the tips, and the tips sometimes additionally rounded; both of these conditions 
have the effect of tending to decrease the contraction of the streams on either 
side of the disk. For circular pipes and disks, the area “a” becomes a 
cresent-shaped opening, and the flow conditions are no longer two-dimension- 
al. Also, the disk angle at closure may be less than 90°; the 4"' and 6" CDC 
valves reported by the authors had closure angles of 75° and 77.59, respec- 
tively. The “model” also assumes no stops, recesses, or projections on the — 
walls, and assumes that the tips of the disk, when at the full closure position, 
contact but are not recessed into, the valve walls. In addition, the two- 
dimensional measurements listed in the authors’ Table I indicate an unequal 
distribution of discharges around the two tips of the disk. : 

Head losses predicted from the “model” would thus not be expected to fur- 
nish very good agreement with actual data at either extreme of disk angle. 
When the valve is fully open, the thickness of the disk and/or disk hub become 
predominate in the formation of the flow geometry; at, or near, full closure 
the calculated area ratio will obviously depart from the actual for valves hav- 
ing maximum @ values of iess than 90°. : 

4 
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In Table I are given the calculated values of Ky, for the “model” over the 
range of disk angles of 00-909. Also given in the table are the values of Ky, 
obtained by the authors for the 4" and 6"' CDC valves, values based upon test 
data by Dickey and Coplen and by Netsch and Schulz (authors’ references (10) 
and (9), respectively), as well as the older values given by Weisbach. The 
values of the discharge coefficient, Ca, as given by the authors, may be con- 
verted to those for Ky by equating comparable expressions for the head loss. 
This yields the relationship, 


. 1.24 
Ky i" pac 


(4) 


he values of Ky for the authors’ references (9) and (10) were recalculated 
(rom the original papers. The values given are not precise, as they were de- 
termined from small-scale test data curves. 

_ As one of the series of check and control valves designed for specialized 
pplications tested recently at the Hydrodynamics Laboratory of the 
ssachusetts Institute of Technology, some characteristics of the butterfly 
ing check valve “H” shown in Fig. 2 were determined. This valve differed 
irom the others considered in this discussion in that there was a smooth 
ransition from an 8" diameter pipe upstream from the valve to an 8" square 
low passage within the valve proper, followed again by a smooth transition to 
p circular pipe. The rectangular disk, allowed to rotate freely, was eccen- 
ically hinged as shown. For steady flows through the valve in the normal, 
f forward, direction, the angles @ assumed by the disk over the range of dis- 


tharges of 1200 to 2000 GPM were 15° and 15930' for two slightly different 
haped disks; the average values of Ky, were 0.82 and 0.66, respectively. For 
ores between 500 and 1200 GPM, the average value of Ky, for the two 
jisks were 0.65 and 0.72, respectively, indicating changes in the internal flow 
4eometry with Reynolds number, these changes in this case possibly being 
‘aused more by the transition section effects upon the flow than by flow 
thanges due to the valve disk. A description of the test facilities and results 


| Some of the data listed in Table I are plotted on Fig. 3. The solid line is 
Yesignated as the “simple model theory.” The test points for swing check 
alve “H” are also plotted. It might be mentioned that data for some com- 
mercial valves of 10" and 24" diameter with circular disks also fall within 
ue spread of the curves plotted. As expected, the differences between the 
imodel” values and the test data are greatest for the largest and smallest 

- The known physical size and shape of the valve and its disk for a particular 
lve could be used to determine more precisely the actual flow area ratios 

nd to give better flow deflection angles in choosing values of C,. It is likely 
at by these more refined calculations the actual Ky, values for valves having 
irelatively simple interior could be estimated more precisely, dependent 
imarily, of course, upon satisfactory values of Cc. The simple “model” 
esented here is not to be considered a rigorous approach to the problem, 

d undoubtedly compensating factors in some cases cause its agreement with 
st data to be fortuitous. The complete characteristics of a given valve could 
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be determined only by a model test. However, the results clearly indicate 
that the butterfly valve is one type of control device where the application of 
the simple energy-momentum relationships will lead to a decent approxima- 
tion of head loss characteristics. The same approach may be used with a 
varying degree of success with other more common types of valves; i.e., gate 
plug, etc. In smaller gate valves especially the head loss coefficient is found 
to vary with the size of the valve because of the relatively greater signifi- 
cance of the boundary configuration, gate recesses, and relatively thicker 
gate disk. Comparisons of various valves, using the area ratio a/A as the 
basis of comparison, are given in (D-3). 
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BUTTERFLY CHECK VALVE ‘H’ 


1348-54 HY 4 August, 1957 


12" Square Disk 


ee 10 Cm. Circ. Disk 


An 
/ 

Simple Model Theory hal 
eae 6" CDC, Circ. Disk =) 


/ 
/ 
I; 
Y), 


Disk Angle « 
FIG. 3 


DISCUSSION 1348-55 


ASCE 


co TT9°0 000°0 906 
00s os77 OOL ‘TT 219°0 ST0°0 008 
OSL BSE OS77 oge O9€ 9L9 819°0 090°0 0k 
Slime 66 16 rs) LY OTT LIT Z€9°0 V/aete) 909 
hd Se ak (4 rad 9 9% 6°tE L°8e 7°0€ LS9°0 EZ" 00S 
g°0OIT Le*6 2°8 T°6 ¢°aT 7°IT 68°6 06°6 SL9°O LGE°O 007 
T6°€ WS°E 9°% 12s L°€ 06°€ 7S°€ Ves 9TL°O 00S°0 oe 
75°T . 78°T 66°0 er e°t Ov°T LTT 76°0 €LL*°O 859°0 90 
Zs°0 8 S7"0 SE°0O YEO Z£°0 7L°0 $9°0 9T°0 z98°0 9z8°0 oot 
00°0 000°T 000°T me) 
xeul Xeul xXeul- 

a=" Wau Rawr ° Sinan Saree ee ee 
qeqnorzg ‘3004 zeqnSaq9 ‘wo OT *ozz9 uc/T €T exenbs yZT 9 9 9@ w7 = ° OTRO D ‘Bay early eATeA 

“TOEGS TOM ah ages Wepdoy pur Kay Ta TTE "7S WOSTSUIOW TSpoR STANTS 


seate, ATZ1a73Ng snoTxte, uO sysez 
IO} 7 eTHuy eaten *sa rly fyuepoTZjze0g ssoy peay 


I AIavl 


vee) 
uh <<. 


; Si) 1} é nie 
; eaeindy 
te acl oy ia¥ of Fe 


. } mate YE 3 ms 
emis het : SEs 


wig), bie ois ee 


4 


. > 
$.i 
vf 
t~ 
an 
r 4 
/ 
as 
fyot 
* eae ' 


ASCE 1348-57 


| Discussion of 
“A STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS” 


by M. B. McPherson and M. L. Karr 
(Proc. Paper 1266) 


CORRECTIONS—As the result of a printer’s error, the illustrations that 
/accompanied Proc. Paper 1266 were the illustrations that were used in Proc. 
Paper 1167. The illustrations for Proc. Paper 1266 are reproduced herewith: 
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SURGE HEIGHTS FOR BUCKET-TYPE ENERGY DISSIPATORS 
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SLOPE OF 1°ON 2 


0.8 | 
A q RANGE OF AS. Vl 
LEGEND 
Yo ni? 


0.3 


0.2 = — é, 4 
SATA RRRERE 


ry ! 
Ay Slew Fai htaes Conese 
Dee eeT naa THiMNAA NAcc ei ATOLL 
Arcee cE TL | | 
te) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
hs | 
hy 


a] 
Cc 
A 
NJ 


PROCEEDINGS PAPERS 


: The technical papers published in the past year are identified by number below. Technical- 
division sponsorship is indicated by an abbreviation at the end of each Paper Number, the 
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to the appropriate issue of “Civil Engineering.” Beginning with Volume 82 (January 1956) papers 
__ were published in Journals of the various Technical Divisions. To locate papers in the Journals, 
the symbols after the paper numbers are followed by a numeral designating the issue of a 
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1061(WW4), 1062(WW4), 1063(WW4), 1064(SU2), 1065(SU2), 1066(SU2)°, 1067(ST5)°, 1068 
(wwa4)°, 1069(Ww4). 


OCTOBER: 1070(EM4), 1071(EM4), 1072(EM4), 1073(EM4), 1074(HW3), 1075(HW3), 1076(HW3), 
| 1077(HY5), 1078(SA5), 1079(SM4), 1080(SM4), 1081(SMm4), 1082(HY5), 1083(SA5), 1084(SA5), 
~ 1085(SA5), 1086(PO5), 1087(SA5), 1088(SA5), 1089(SA5), 1090(HW3), 1091(EM4)°, 1092 
| _ (Hy5)°, 1093(HW3)°, 1094(Po5)¢, 1095 (Sma4)°. 


NOVEMBER: 1096(ST6), 1097(ST6), 1098(ST6), 1099(ST6), 1100(ST6), HOSTS), 1102(IR3), 1103 
(IR3), 1104(1R3), OS(IR3), 106(ST6), 1107(ST6), LOS(ST6), O9(AT3), W10(AT3)°, 11(IR3)°; 
_ miasre)¢. 


DECEMBER: 1113(HY¥6), 1114(HY6), 1115(SA6), 1116(SA6), 1117(SU3), 1118(SU3), 1119(WW5), 
+ 1120(WW5), 1121(WW5), 1122(WW5), 1123(WW5), 1124(WW5)°, 1125(BD1)°, 1126(SA8), 1127 
 (SA6), 1128(WW5), 1129(SA6)°, 1130(PO6)°, 1131(H¥6)°, 1132(PO6), 1133(PO6), 1134(PO6), 
_ 1135(BD1). 


VOLUME 83 (1957) 


| ARY: 1136(CP1), 1137(CP1), 1138(EM1), 1139(EM1), 1140(EM1), 1141(EM1), 1142(SM1), 
| 1143(SM1), 1144(SM1), 1145(SM1), 1146(ST1), 1147(ST1), 1148(ST1), 1149(ST1), 1150(ST1), 

 -1151(ST1), 1152(CP1)©, 1153(HW1), 1154(EM1)°, 1155(SM1)°, 1156(ST1)°, 1157(EM1), 1158 
| _ (®M1), 1159(SM1), 1160(SM1), 1161(SM1). 


BRUARY: 1162(HY1), 1163(HY1), 1164(HY1), 1165(HY1), 1166(HY1), 1167(HY1), 1168(SA1), 

ii oo : cement reakie 1172(SA1), 1173(SA1), 1174(SA1), 1175(SA1), 1176(SA1), 
arene, 1178(SA1), 1179(SA1), 1180(SA1), 1181(SA1), 1182(PO1), 1183(PO1), 1184(PO1), 
1185(PO1)°. 


MARCH: 1186(ST2), 1187(ST2), 1188(ST2), 1189(ST2), 1190(ST2), 1191(ST2), 1192(ST2)°, 1193 
 (PL1), 1194(PL1), 1195(PL1), A ee Sad 
2IL: 1197(HY2), 1198(HY2), 1199(HY2), 1200(HY2), : HY2), 1 
(SA2), neuen, peels Pouens. 1207(SM2), 1208(WW1), 1209(WW1), 1210(ww1), 
1211(WW1), 1212(EM2), 1213(EM2), 1214(EM2), 1215(PO2), 1216(PO2), 1217(PO2), 1218 
 (GA2), 1219(SA2), 1220(SA2), 1221(SA2), 1222(SA2), 1223(SA2), 1224(SA2), 1225(PO)°, 1226 
- (Ww1)¢, 1227(SA2)°, 1228(SM2)°, 1229(EM2)°, ee EA, sere wuntauaeten 
AY: , 1232(ST3), 1233(ST3), 1234(ST3), 1235(IR1), 12: ‘ , 1238(WW2), 
Saris. eunvial: Bory aig 1242(WW2), 1243(WW2), 1244(HW2), 1245(HW2), 1 
|| (Hw2), 1247(HW2), 1248(ww2), 1249(HW2), 1250(AW2), 1251(WW2), 1252(WWw2), 1253(0R1), 
_ 1254(S73), 1255(ST3), 1256(HW2), 1257(IR1)°, 1258(HW2)°, 1259(ST3)°. | 
‘UNE: 126 ; , 1262(H¥3), 1263(HY3), 1264(HY3), 1265(HY3), 1266(HY3), 1267 
; se ete ate mia), 1271(SA3), 1272(SA3), 1273(SA3), 1274(SA3), 1275 
_ (SAS), 1276(SA3), 12 77(HY3), 1278(HY3), 1279(PL2), 1280(PL2), 1261(PL2), 1282(SA3), 1283 
_ (HY3)%1284(PO3), 1285(PO3), 1286(PO3), 1287(PO3)”, 1288 a pel ae 
ULY:. 1290(EM3) , 1292(EM3), 129 , 1294( » 1295 ; 
ate = os Se ee », 1301(SM3) , 1902(ST4) ,_ 130; 
- (ST4), 1304(ST4), 4), 1909(SM3), 1310(SUt) ©, 
--1311(EM3)°, 1 1314(ST4), 1315(ST4) , S74), 1317(ST4), 1318 
ST4), 1319(SM3)° 1920(ST4) , 4), 1922(EM3) , 1925(AT1), 1924(AT1), 1925(ATI), 
33(SA4), 1334(SA4), 1335(SA4), 1336(SA4), 
"sH(co!), 142(CO4}, 1343(CO1), 1344(PO4), 
, 1349(SA4)", 1350(PO4), 1351(PO4), 
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